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SUMMARY 


Metastable  molecular  fuels,  high-energy  long-lived  molecular  species  that  do  not 
decay  by  radiation,  tunneling,  or  other  means  when  isolated  in  vacuum,  have  been 
proposed  as  the  basis  for  possible  new  propulsion  schemes.  This  theoretical  project  has 
been  devoted  to  examination  of  ion-pair  species  bound  by  the  coulomb  attraction  between  a 
stable  positive  and  stable  negative  ion  as  possible  candidate  fuels.  Theoretical  results  for 
H3O,  U3H,  and  H3F  are  presented  in  this  report  along  with  specific  impulse  predictions 
based  on  calculated  energies. 

The  greatest  effort  was  directed  at  the  H3O  molecule,  which  was  selected  as  the  first 
candidate  ion-pair  system  because  of  the  significant  electron  affinity  (1.46  eV)  of  the 
oxygen  atom.  Minimum  geometries  and  transition  states  on  the  first  two  doublet  surfaces 
of  H3O  have  been  obtained  by  multiconfiguration  self-consistent  field  (MCSCF)  analytic 
gradient  techniques,  and  the  correlation  diagram  connecting  these  states  to  the  separated 
fragment  asymptotes  to  which  they  fragment  has  been  constructed.  While  the  energetic 
ion-pair  state  was  found  unstable,  the  possible  stability  of  the  lowest  state  of  H3O  has  been 
the  subject  of  significant  experimental  interest.  Theoretical  results  have  been  related  to 
experimental  reports  of  metastable  D3O. 

While  there  exists  a  region  of  the  lowest  energy  surface  that  may  be  described  as  an 
ion-pair,  our  results  have  shown  that  the  ion-pair  is  not  even  a  local  minimum  on  that 
surface.  The  ion-pair  minimum  identified  as  a  local  minimum  on  the  excited  1  2E  surface  in 
C3V  symmetry  is  not  stable  in  lower  symmetry  and  is  predicted  to  dissociate  without  an 
energy  barrier  to  H2  +  H  +  0.  This  observation  has  been  qualitatively  explained  by  back¬ 
charge  transfer  from  O*  to  H3 ,  leading  to  neutral  H3,  which  Is  unstable  with  respect  to 
H2  +  H. 

In  contrast,  the  Li3H  molecule  is  more  likely  to  exist  as  a  stable  ion-pair  because 
both  Li|  and  U3  are  bound  species.  As  anticipated  from  this  argument,  the  ion-pair 
conformation  is  found  to  be  a  local  minimum  with  all  real  frequencies  on  the  lowest 
potential  surface  of  U3H.  However,  there  exists  a  global  planar  minimum  20  kcal/mol 
lower  in  energy  and  the  barrier  height  for  conversion  of  the  ion-pair  state  to  the  planar  form 
is  predicted  to  be  only  1.3  kcal/mol.The  correlation  diagram  for  Li3H  has  been  determined. 
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In  contrast  to  H3O,  the  lowest  states  of  L13H  lie  below  all  possible  fragment  energies. 

While  the  lowest  state  of  LyH  is  most  definitely  stable,  it  is  not  particularly  energetic. 

Calculations  of  the  minima  and  frequencies  of  the  ion-pair  state  of  H3F  support  the 
conclusion  that  ion-pair  states  based  on  the  H3  cation  are  not  stable  with  respect  to  back¬ 
charge  transfer. 

The  goal  of  this  project  was  to  identify  energetic  molecular  states  that  might  serve  as 
the  basis  of  new  propulsion  schemes.  To  estimate  the  performance  of  these  new  metastable 
molecules,  the  specific  impulse  (Isp)  was  calculated  for  each  metastable  molecule.  The  1^ 
calculations  assumed  that  metastable  molecules  (1)  could  be  stabilized,  (2)  used  the  same 
exhaust  conditions,  and  (3)  did  not  optimize  Isp.  Metastable  D3O  as  a  monopropellant  is 
predicted  to  provide  a  7%  improvement  over  H2/O2  under  the  same  conditions.  Because  the 
ground  state  of  LyH  lies  more  than  30  kcal/mol  below  LiH  +  U2,  LvjH  as  a  fuel  does  not 
have  a  favorable  1^,.  However,  U3H  as  an  additive  to  H2  in  the  ratio  Li3H:H2:02  of 
5:30:4  is  predicted  to  provide  a  16%  improvement.  Use  of  L12H  as  an  additive  is  estimated 
to  produce  a  20%  improvement. 
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INTRODUCTION 


The  energy  density  of  today's  rocket  propellants  severely  limit  rocket  payloads  and 

*  mission  capability.  For  example,  no  matter  how  efficient  the  engine,  a  single  stage  rocket 
will  always  have  a  negative  payload-  Therefore,  novel  concepts  that  might  serve  as  the 
basis  of  new  propulsion  schemes  are  of  great  interest  One  concept  that  has  recendy  been 
considered  is  the  use  of  metastable  molecular  fuels.  Metastable  molecules  are  long-lived 
molecular  species  that  do  not  decay  by  radiation,  tunneling,  or  other  means  when  isolated 
in  vacuum. 

The  specific  impulse,  Isp,  which  is  proportional  to  the  square  root  of  the  ratio  of  the 
heat  of  reaction  to  the  molecular  weight  of  the  products,  is  the  figure  of  merit  for  evaluating 
propulsion  reactions.  Therefore,  ideal  rocket  propellants  are  compounds  with  high  energy 
content  which  form  low  molecular  weight  products.  This  requirement  virtually  restricts  the 
candidate  species  for  metastable  molecular  fuels  to  molecules  composed  of  first-row  atoms. 

Theoretical  prediction,1  several  years  ago,  of  an  extremely  energetic  state  of  H4 
initially  excited  great  interest  With  very  high  energy  and  low  molecular  weight,  H4 
appeared  to  be  a  promising  molecular  fuel.  However,  detailed  theoretical  study  by  several 
groups24  during  the  last  two  years  has  led  to  the  conclusion  that  the  predicted  state  of  H4  is 
likely  to  have  an  extremely  short  lifetime.  Nonetheless,  the  serious  examination  of  the  H4 
system  has  provided  valuable  insight  that  can  guide  the  selection  and  study  of  future 
candidate  molecules. 

Tetrahydrogen  is  only  one  example  of  a  class  of  excited  compounds  that  might  be 
termed  ion-pair  states  because  the  stable  geometry  is  predicted  to  resemble  an  H3  ion  in  an 
equilateral  triangle  with  the  H"  ion  above  the  center.  The  fate  of  an  ion-pair  species,  once 
formed,  will  depend  sensitively  on  the  details  of  the  potential  surfaces  that  govern  decay 
processes  such  as  optical  transitions,  predissociation,  and  internal  conversion.  Thus,  it  is 
worthwhile  to  investigate  various  potential  ion-pair  energetic  species  to  find  the  case  that  is 
most  stable.  This  report  presents  results  of  theoretical  studies  of  candidate  ion-pair  species 
composed  of  first-row  atoms. 

* 

We  expect  the  strength  of  the  ion-pair  bond  to  depend  on  the  ability  of  the  anion  to 
remain  negatively  charged,  i.e.,  on  the  electron  affinity  of  the  corresponding  neutral. 
Because  the  oxygen  atom  has  an  electron  affinity5  of  1.46  eV,  nearly  twice  that  of  atomic 
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hydrogen,  H3O  was  selected  as  the  first  candidate  system.  Significant  effort  was  devoted  to 
this  system,  with  the  two-fold  goal  of  (1)  characterizing  as  completely  as  possible  the 
minima  and  barriers  on  the  low-lying  potential  surfaces  of  this  prototype  ion-pair  system 
and  (2)  developing  and  validating  theoretical  strategies  that  can  be  applied  to  other  systems. 
While  the  ion-pair  state  of  H3O  is  not  predicted  to  be  stable,  our  study  examined  other 
portions  of  the  potential  surface,  and  addressed  the  metastability  of  the  lowest  state 
(12A')-6*11 

Two  other  ion-pair  state  molecules,  U3H  and  H3F,  were  also  examined.  We  have 
considered  the  U3H  system  as  an  example  of  an  ion-pair  based  on  a  different  cation.  In 
addition  to  addressing  the  general  question  of  when  can  we  expect  the  coulomb  interaction 
to  dominate  molecular  bonding,  our  work  on  U3H  has  proved  relevant  to  the  suggestion 
that  addition  of  Li  would  improve  the  hydrogen-oxygen  fuel  if  the  technical  problem  of  Li 
metal  formation  could  be  overcome.  Since  fluorine  has  the  largest  electron  affinity  of  any 
first  row  element  (3.4  eV),5  the  H3F  ion-pair  state  was  also  evaluated. 

This  report  summarizes  the  results  of  ab  initio  calculations  on  H3O,  L^H,  and  H3F 
and  presents  some  general  conclusions  with  respect  to  stability  of  ion-pair  species.  Details 
of  the  theoretical  calculations  are  given  in  the  resulting  publications,  which  are  attached  as 
appendices.  The  specific  impulse  to  be  expected  for  H3O  and  L^H,  the  former  as  a 
monopropellant  and  the  latter  as  an  additive  to  H2,  has  been  predicted  based  on  the 
calculated  heats  of  formation.  Tables  which  summarize  the  theoretical  performance  of  these 
fuels  are  provided  in  the  section  entitled  Calculations  of  Specific  Impulse. 
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CALCULATIONS  OF  ION-PAIR  STATES 


This  theoretical  project  was  devoted  to  study  of  ion-pair  species  bound  by  the 
coulomb  attractic'  between  a  positive  and  negative  ion  as  possible  candidate  fuels.  The 
greatest  effort  was  directed  at  the  H3O  molecule,  which  was  selected  as  the  first  candidate 
system  on  the  basis  of  the  significant,  electron  affinity5  (1.46  eV)  of  the  oxygen  atom. 
Results  obtained  for  H3O,  U3H,  and  H3F  systems  are  summarized  below. 

POTENTIAL  SURFACES  OF  H3O 

Features  of  the  two  lowest  doublet  potential  surfaces  of  H3O,  minima  and  transition 
states,  are  illustrated  in  Figure  1,  in  which  the  separated  fragment  asymptotes  to  which  they 
dissociate  are  also  indicated.  Geometries  have  been  optimized  by  multiconfiguration  self- 
consistent  field  (MC3CF)  analytic  gradient  techniques,  and  large-scale  multireference 
configuration  interaction  (Cl)  calculations  have  been  used  to  obtain  energy  separations. 

The  calculations  and  conclusions  are  described  in  detail  in  a  manuscript12  that  has  been 
published  in  the  Journal  of  Chemical  Physics  and  is  attached  as  Appendix  A. 

While  there  exists  a  region  of  the  lowest  energy  surface  that  may  be  described  as  an 
ion-pair,  our  results  have  shown  that  the  ion-pair  is  not  even  a  local  minimum  on  that 
surface.  The  ion-pair  minimum  identified  on  the  excited  1%  surface  in  C^  geometry 
(shown  by  a  dashed  line  in  Figure  1)  is  a  local  minimum.  It  is,  however,  not  stable  in 
lower  symmetry  and  is  predicted  to  dissociate  without  an  energy  barrier  to  H2  +  H  +  O. 
This  observation  has  been  qualitatively  explained  by  noting  that  the  neutral  H3  is  not  a 
stable  species.  Therefore,  any  back-charge  transfer  from  O'  to  H3  will  lead  to  the  neutral 
H3,  which  is  unstable  with  respect  to  H2  +  H. 

In  agreement  with  previous  work,1 3  we  find  the  lowest  state  of  H3O  may  be 
characterized  as  an  HjO4-  core  surrounded  by  an  oxygen  3s  Rydberg  electron.  This  finding 
supports  the  inclusion  of  the  H3O  molecule  in  a  recent  discussion  by  Herzberg14  of 
hypervalent  hydrides  that  may  be  expected  to  exist  as  Rydberg  molecules  because  of  the 
stability  of  the  corresponding  ion.  As  shown  in  Figure  1,  the  energy  of  the  lowest  state  of 


5 


ENERGY  (KcaV  mol) 


Figure  1.  Minima  and  transition  states  for  the  H3O  system  obtained  in  this  work. 

Energies  are  given  relative  to  H2O  +  H.  The  12E  ion-pair  state  (dashed) 
dissociates  without  a  barrier  to  H2+H+O.  The  dashed  transition  state  on  the 
1  2A"  potential  (labeled  valence)  surface  connects  to  the  OH+H+H 
asymptote  (not  shown). 
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H3O  is  clearly  above  that  of  the  asymptotes,  OH  +  H2  and  H2O  +  H.  The  stability  of  the 
H30(12A')  state  depends  on  the  barriers  in  the  potential  energy  surface  to  dissociation. 

The  present  work  confirms  the  understanding  obtained  previously,12  i.e.,  that  the 
vibrationless  potential  of  HsO^A')  has  a  very  small  barrier,  3.58  kcal/mol  according  to 
the  present  calculations,  for  dissociation  to  H2O  +  H.  Including  zero-point  energy,  a 
barrier  of  0.4  kcal/mol  is  predicted  for  this  dissociation.  Using  the  ratio  of  masses  to 
estimate  the  vibrational  frequencies  for  D3O,  we  estimate  a  barrier  of  1.3  kcal/mol  for 
dissociation  of  D3O. 

There  has  been  significant  experimental  interest6-11  in  characterizing  the  H3O 
system  and  in  establishing  the  metastability  of  the  lowest  state.  Although  an  earlier  mass 
spectrometric  investigation^  reported  a  lifetime  of  at  least  1  p.s  for  both  H3O  and  D3O,  the 
work  of  March  and  Young,9  by  a  very  similar  experimental  technique  leads  to  the 
conclusion  that  H3O  is  not  stable.  Gellene  and  Porter6  have  reported  observation  of 
metastable  D3O  with  a  lifetime  >  0.6  its  but  found  no  evidence  of  H3O.  Our  prediction  of  a 
1.3  kcal/mol  dissociation  barrier  for  D3O,  but  only  a  0.4  kcal/mol  barrier  for  H3O,  may  be 
consistent  with  this  report.  The  energy  of  the  lowest  state  of  H3O  with  respect  to  H2O  + 

H  is  being  probed  experimentally  by  Devynck  and  Peterson10  in  this  laboratory.  (These 
workers  are  also  investigating  the  D3O  isotope.)  Preliminary  results  are  in  good  agreement 
with  the  23  kcal/mol  excitation  energy  of  H3O  above  H2O  +  H  predicted  in  this  work.  A 
more  complete  understanding  of  this  system  is  expected  from  an  experiment  with 
photoionization  detection  which  will  establish  the  energy,  with  respect  to  the  ion,  of  the 
state  being  formed.  Such  experiments  are  also  under  way  in  this  laboratory.10-11 

CALCULATIONS  ON  L13H 

As  discussed  specifically  for  H3O,  back-charge  transfer  is  now  thought  to  limit  the 
stability  of  any  ion-pair  species  with  H3  as  the  positive  ion.  Michels15  has  suggested  that 
Li3H  may  be  a  more  attractive  candidate  because  both  Li|  and  Li3  are  stable  species. 
Preliminary  calculations  by  Michels  and  Montgomery15  support  the  stability  of  ground  state 
Li3H.  We  therefore  studied  the  low-lying  states  of  LyH  using  a  similar  computational 
procedure  to  that  followed  for  H3O  and  established  the  correlation  diagram  linking  L^H 
states  to  the  LiH  +  Li2,  Li2H  +  Li,  and  LI3  +  H  asymptotes. 
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In  agreement  with  a  preliminary  report15  and  with  the  recent  study  of  Montgomery 
et  al.,16  we  find  a  local  minimum  with  all  real  vibrational  frequencies  on  the  lowest 
potential  surface,  the  l^A’,  at  the  Csv  symmetry  pyramidal  geometry.  Details  are  provided 
in  the  resulting  publication  attached  as  Appendix  B.17  Analysis  of  the  wavefunction 
indicates  substantial  mixing  of  ion-pair  and  valence  character  at  this  pyramidal  geometry. 
The  Li-Li  bond  length  of  5.173  ao  for  U3H  is  ~0.7  ao  shorter  than  the  corresponding  bond 
length  in  IJ3,  in  contrast  to  the  case  for  H3O,  where  the  ion-pair  state  was  found  at  the  H- 
H  separation  of  H|.  The  Li-Li  bond  lengths  in  the  C2v  structure  of  Li3  (2Ai),  to  which 
this  state  correlates,  are  longer  as  well 

More  important,  however,  there  exists  a  planar  minimum  20.3  kcal/mol  lower  in 
energy  on  the  lowest  singlet  potential  surface,  as  shown  in  Figure  2.  We  also  optimized 
the  transition  state  for  conversion  of  the  pyramidal  form  to  the  planar  geometry  that  is 
predicted  to  have  a  low  imaginary  frequency  of  120  i  cm-1.  The  transition  occurs  at  a  very 
modest  geometrical  distortion  from  the  pyramid,  with  a  calculated  barrier  height  of  1.3 
kcal/mol  including  zero-point  energy.  This  result  is  in  agreement  with  the  recent  work  of 
Montgomery  et  al.,16  who  found  the  transition  state  at  values  between  <1  kcai/mol  above 
and  0.05  kcal/mol  below  the  pyramidal  minimum.  These  results  include  zero-point  energy , 
and  depend  on  the  basis  set  and  level  of  calculation). 

The  stability  of  the  L^H  (UA')  ion-pair  should  be  considered  to  be  problematical 
because  of  this  low  or  nonexistent  barrier.  The  lowest  triplet  state  has  a  planar  minimum 
geometry  that  lies  17.9  kcal/mol  above  the  singlet  ground  state.  The  triplet  does  not  have  a 
local  minimum  in  a  pyramidal  geometry. 

As  shown  in  Figure  2,  in  contrast  to  the  case  for  H3O,  the  singlet  and  triplet  states 
of  Li3H  lie  below  all  possible  separated  fragment  asymptotes.  Figure  2  also  includes  the 
second  state  of  each  symmetry  at  the  geometry  of  the  ground  state.  The  energies  are  of 
qualitative  reliability  only.  These  states  are  predicted,  however,  to  be  bound  with  respect  to 
the  LiH  +  Li2  asymptotes  to  which  they  correlate. 

The  prediction  of  a  true  minimum  on  the  12A'  surface  (all  real  frequencies)  at  the 
Qjv  geometry  of  mixed  ionic  and  covalent  character  supports  the  original  suggestion  that 
because  of  the  stability  of  Li3  and  Li3,  back-charge  transfer  will  not  limit  the  stability  of 
LyH  with  respect  to  dissociation.  Unlike  the  H4  and  H3O  cases,  there  is  in  this  molecule  a 
lower  energy  planar  form.  Given  the  very  low  calculated  barrier  to  interconversion,  the 
pyramidal  local  minimum  is  predicted  to  be  unstable  relative  to  the  global  planar  minimum. 
While  the  lowest  state  of  Li3H  is  most  definitely  stable,  it  is  not  particularly  energetic. 
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LiH  +  LI2  Pyramidal 


Planar  Ll2H  +  LI 


UH  +  L^lhlg) 
UH  +  U2(B’nu} 
uhTu^zJ 


i'a- 

Li3(2A,)  +  H(2S)  /  i3a. 


Li2H(2r)  +  U(2S) 


Correlation  diagram  linking  the  states  of  U3H  to  those  of  Li3  +  H  and 
U2H  +  H,  drawn  in  Cs  symmetry  using  calculated  energies.  The  second 
state  of  each  symmetry  has  been  calculated  at  the  geometry  of  the  lower 


ION-PAIR  STATE  OF  H3F 


The  goal  of  the  calculations  on  H3F  was  to  conclusively  demonstrate  whether  back¬ 
charge  transfer  forming  the  dissociative  ground  state  of  H3  will  always  prevent  stability  of 
an  ion-pair  state  with  H3  as  the  cation,  as  indicated  by  the  information  on  H3O  and  H4. 
Since  fluorine  has  the  largest  electron  affinity  of  any  first-row  atom  (3.4  eV)5,  fluorine  was 
expected  to  be  the  most  favorable  anion.  Calculations  were  limited  to  determining  the 
optimized  geometry  and  harmonic  frequencies  of  the  ion-pair  state  of  H3F  and  a  brief 
search  for  the  global  minimum  on  the  lowest  potential  surface.  In  order  to  obtain  results  as 
rapidly  as  possible,  preliminary  optimization  was  performed  at  the  self-consistent  field 
(SCF)  level. 

Optimized  geometries  and  frequencies  are  summarized  in  Table  1.  The  F  atom  basis 
set  was  taken  from  the  work  of  Bauschlicher  and  Taylor1**  on  the  electron  affinity  of  the 
fluorine  atom.  The  qualitative  conclusions  presented  here  are  independent  of  the  basis  set 
and  type  of  calculation.  The  C3V  symmetry  pyramidal  minimum  of  the  ion-pair  state, 
which  is  found  for  an  H-H  separation  very  close  to  that  of  H3  has  a  doubly  degenerate  E- 
symmetry  imaginary  frequency.  At  the  SCF  level,  distortion  of  this  geometry  according  to 
the  normal  modes  corresponding  to  these  imaginary  frequencies  leads  to  optimization  of  a 
very  slightly  bound  van  der  Waals  complex  of  H2  and  HF.  Two  minima  very  close  in 
energy  are  found;  in  one  of  them  the  H2  is  attached  to  the  hydrogen  end  of  HF  and  in  the 
other  to  the  fluorine.  The  H2--FH  result  is  in  reasonable  agreement  with  the  theoretical 
result  reported  by  S apse. 19  The  van  der  Waals  complexes  are  predicted  to  be  bound  by  <1 
kcal/mol  at  the  SCF  level.  They  have  no  energy  content. 

Although  the  study  was  not  carried  through  to  convergence  at  the  MCSCF  level, 
similar  results  are  indicated.  As  shown  in  Table  1,  at  the  MCSCF  level  the  C3v  symmetry 
ion-pair  state  has  a  doubly-degenerate  imaginary  frequency  in  agreement  with  the  SCF 
prediction.  Preliminary  results  indicate  that  MCSCF  optimization  will  lead  to  the  van  der 
Waals  complex,  in  agreement  with  ihe  SCF  result.  No  stable  ion-pair  state  of  H3F  could 
be  found. 

The  calculations  on  H3F  reported  here  support  the  conclusion  that  the  H3  cation 
does  not  form  stable  ion-pair  states. 
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Table  1. 
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CALCULATIONS  OF  SPECIFIC  IMPULSE 


The  goal  of  this  project  is  to  identify  energetic  molecular  states  that  might  serve  as 
the  basis  of  new  propulsion  schemes.  It  is  therefore  valuable  to  estimate  the  specific 
impulse  (Isp)  that  could  be  expected  from  the  species  investigated  here,  under  the 
assumption  the  metastable  states  can  be  stabilized.  For  these  estimates,  the  NASA-Lewis 
code2*)  for  prediction  of  rocket  performance  has  been  used  to  calculate  specific  impulse 
without  extensive  optimization. 

H30 

Results  for  various  ratios  of  H3O  or  D3O  and  O2  as  well  as  for  H2  +  1/2  O2  for 
comparison  are  shown  in  Table  2.  A  heat  of  formation  of  12.458  kcal/mol  for  the  lowest 
state  of  H3O  (or  D3O)  was  used  in  these  calculations.  This  value  was  computed  from  our 
calculated  relative  energy  with  respect  to  H2O  +  H  and  literature  values21  for  the  heat  of 
formation  of  H2O  and  H.  The  unstable  energetic  ion-pair  state  was  not  considered.  Under 
the  conditions  considered  thus  far,  H2  +  1/2  O2  is  predicted  to  have  an  Ijp  of  369  sec, 
whereas  values  greater  than  400  sec  are  often  quoted.  D3O,  as  a  monopropellant,  is 
predicted  to  be  more  favorable,  with  an  1^  of  396  sec.  The  lighter  H3O  (which  is, 
however,  not  predicted  to  be  metastable),  would  still  be  more  favorable,  with  a  predicted 
value  of  427  sec. 
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Table  2. 


ESTIMATED  1^  FOR  H30  AND  D3O 


Isp(s) 

h2  +  j02  369.3 

H3O+I02  351.9 

H30+j02  372.9 

H3O  427.0 

D30+|o2  337.3 

D30+io2  354.9 

°30  395.7 
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Li3H  AND  RELATED  COMPOUNDS 


Using  the  calculated  electronic  energies,  zero  point  energies,  and  literature  values 
for  the  heats  of  formation2*  of  LiH  and  Li2,  we  calculate  the  heats  of  formation  of  the  Li3H 
*A'  planar  global  minimum  and  ion-pair  (pyramidal)  configurations  to  be  51.06  and  71.36 
kcal/mol,  respectively.  The  heat  of  formation  of  the  3A"  state  is  predicted  to  be  69.62 
kcal/mol  In  Table  3,  we  list  the  specific  impulse  calculated  with  the  NASA-Lewis  code20 
for  heats  of  formation  approximating  these  values  for  combinations  of  Li3H,  02.  and  H2. 
The  H2  was  included  following  the  suggestion  of  Rodgers,22  who  pointed  out  the 
improvement  in  predicted  performance  of  H2-O2  with  addition  of  Li.  Ratios  have  been 
only  approximately  optimized. 

Because  the  ground  state  of  LrjH  lies  more  than  30  kcal/mol  below  LiH  +  Li2,  the 
molecule  itself  as  a  fuel,  even  if  the  triplet  state  or  pyramidal  form  of  the  singlet  can  be 
used,  is  predicted  to  provide  a  low  Ijp.  However,  Li3H  looks  promising  as  an  additive  to 
H2,  with  predicted  Isp  values  of  429  and  438  sec,  respectively,  for  assumed  heats  of 
formation  of  68  and  50  kcal/mol  in  the  molar  ratio  Li3H:H2:C>2  of  5:30:4.  Recall  that  under 
these  conditions  H2  +  1/2  O2  has  an  1^  of  369  s.  Finally,  from  our  calculated  results  we 
roughly  estimate  the  heat  of  formation  of  Li2H  to  be  50  kcal/mol.  Use  of  LnH  as  an 
additive  is  estimated  to  be  even  more  favorable,  with  a  calculated  specific  impulse  of  445 
sec. 
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Table  3 

ESTIMATED  %  FOR  Li  COMPOUNDS 
Molar  Composition  _ _ ISp  (s) 


0 1 

2 

1 

369 

LhH 

02. 

AHf  =  68a 

AHf  =  503 

5 

2 

352 

331 

5 

20 

2 

416 

5 

30 

4 

438 

429 

LipH 

H?. 

_Qz 

AHf  =  503 

5 

2 

361 

5 

30 

4 

445 

^  Energies  in  kcal/mol. 
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CONCLUSIONS  AND  IMPLICATIONS  FOR  FUTURE  WORK 


Several  conclusions  may  be  drawn  from  these  studies  of  ion-pair  states.  First,  in 
agreement  with  previous  indications,22  addition  of  compounds  of  Li  is  predicted  to  enhance 
performance  of  the  H2-O2  system.  Entrainment  of  atomic  Li  in  hydrogen  is  made  difficult 
by  the  tendency  of  Li  to  agglomerate  to  form  the  metal.22*23  Introduction  of  Li  in  the  form 
of  LyH  or  Li2H  may  ameliorate  this  problem.  However,  compounds  of  Li  may  also  tend 
to  form  clusters;  formation  of  oligomers  of  LiF,24  for  example,  has  been  the  subject  of 
study.  Thus,  theoretical  studies  of  dimer  formation  of  these  compounds  and  experimental 
investigation  of  stabilization,  most  probably  in  low-temperature  matrices,  appear  to  be 
valuable. 

Second,  from  our  work  on  H3O  and  H3F  and  the  work  of  Michels2*15  and  Lester,3 
among  others,  on  H4  and  H3Li,  the  conclusion  that  charge  separated  species  with  H3  as  the 
positive  ion  are  destabilized  by  back-charge  transfer  seems  to  be  fairly  well  established.  In 
proposing  such  species,  the  fact  that  covalent  bonding  often  competes  energetically  with  the 
ionic  configuration  should  not  be  overlooked.  This  seems  to  be  particularly  troublesome 
for  an  open  shell  case,  such  as  H3O,  in  which  a  number  of  states  arise  from  the  same 
limits,  but  even  H3F,  which  has  only  one  state  arising  from  the  ion-pair  limit  and  no 
significant  covalent  bonding,  does  not  have  a  stable  ion-pair  state.  No  further  effort  seems 
warranted  on  systems  which  include  the  H3  moiety.  However,  it  may  be  possible  to  exploit 
the  large  electron  affinity  of  atomic  fluorine  and  certain  fluorine-containing  compounds, 
(e.g.,  BeF3,  BF4)  by  combining  them  with  cations  for  which  the  corresponding  neutrals 
have  particularly  small  ionization  potentials.  This  would  minimize  the  back-charge  transfer 
mechanism,  which  has  limited  the  stability  of  many  of  the  ion-pair  systems  studied  to  date. 
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The  first  two  doublet  and  quartet  states  of  H3  O  have  been  surveyed  by  multiconfiguration  self- 
consistent  field/first-order  configuration  interaction  (MCSCF/FOCI)  calculations  in  C3„ 
symmetry.  Geometries  of  the  minima  on  the  doublet  surfaces  have  been  optimized  by  MCSCF 
gradient  techniques  and  energies  obtained  by  large-scale  multireference  single  and  double 
excitation  Cl  calculations.  The  correlation  diagram  linking  the  minima  to  different  dissociation 
limits  has  been  established.  A  local  tr-nimum  in  C, ,  symmetry  of  ion-pair  character  is  shown 
to  be  unstable  with  respect  to  dissociation  to  H2  +  H  +  O.  The  lowest  state,  1 2A,  in  C3.. 

( 1 2A '  in  C, )  may  be  characterized  as  an  H3  O  +  core  surrounded  by  an  oxygen  3s  Rydberg 
electron.  Transition  states  for  dissociation  of  the  1  2 A '  state  to  H20  -f  H  and  to  OH  +  H2  have 
been  investigated.  An  extremely  low  barrier  height,  3.58  keal/mol  without  vibration,  0.4 
keal/mol  for  HjO,  an  estimated  1.3  keal/mol  for  D3  0,  with  zero-point  energy,  is  found  for 
dissociation  of  the  1  2A  '  state  to  H20  +  H.  Within  the  uncertainty  of  the  calculation  it  is  not 
possible  to  predict  whether  the  lowest  state  of  H3  O  should  be  observable  experimentally.  The 
relationship  of  these  results  to  experimental  observations  for  H30  is  discussed. 


INTRODUCTION 

There  is  considerable  intrinsic  interest  in  understanding 
novel  bonding  situations  in  polyatomic  molecules.  In  addi¬ 
tion,  stable  states  of  molecules  with  high  energy  content  are 
of  potential  practical  significance  for  the  development  of 
new  propellants.1  A  recent  model  based  on  charge  separa¬ 
tion  in  an  excited  species  has  been  proposed  to  predict  the 
existence  of  stable  excited  states  even  in  compounds  with 
unstable  ground  states  and  has  been  applied  to  explain*  a 
calculated  local  minimum  in  excited  H«.  At  the  geometry  of 
the  minimum,  the  system  may  be  described  as  H3+  in  its 
equilateral  triangle  ground  state  with  a  H"  ion  above  the 
center  of  the  triangle  ar.d,  therefore,  may  be  considered  an 
ion-pair  state.  While  subsequent  calculations3,4  have  indicat¬ 
ed  the  lowest  state  of  H«  can  be  expected  to  decay  to  2  H2, 
the  class  of  ion-pair  states  of  polyatomics  remains  of  general 
interest  and  largely  unexplored. 

We  have  undertaken,  therefore,  a  series  of  theoretical 
surveys  of  ion-pair  and  nearby  states  of  small  polyatomic 
molecules.  In  this  paper,  we  report  on  our  initial  work  on  the 
H30  system.  Because  the  strength  of  the  ion-pair  bond  is 
expected  to  be  in  some  sense  proportional  to  the  stability  of 
the  negative  ion  and  the  oxygen  atom  has  an  electron  affinity 
of  1.46  eV,5  nearly  twice  that  of  atomic  hydrogen,  H30  was 
selected  as  an  interesting  candidate  system.  Furthermore, 
there  is  substantial  theoretical*  and  experimental7*10  interest 
in  the  existence  of  stable  or  metastable  H3  O,  in  its  spectros¬ 
copy1 1,11  as  an  example  of  a  Rydberg  molecule,  and  in  the 
potential  surface13"14  for  the  hydrogen  abstraction  reaction 
0H  +  H2-H20  +  H. 

As  reviewed  in  previous  work,4,7,9  the  possible  existence 
of  the  HjO  radical  has  been  the  subject  of  investigation  for 
more  than  20  years.  Experimental  evidence  for  metastable 
states  of  D3  0,  but  not  of  H3  O,  has  been  reported  by  Gcllene 
and  Porter7  from  analysis  of  fast  neutral  beam  scattering 
profiles  and  collisional  reionization  mass  spectra.  A  subse¬ 


quent  paper  by  Raskit  and  Porter,*  however,  found  no  evi¬ 
dence  of  metastable  D3  ‘*0.  Griffiths  el  al .*  have  reported 
experimental  evidence  by  neutralization-reionization  spec¬ 
troscopy  of  HjO  radicals  with  a  lifetime  of  0.41  fts.  How¬ 
ever,  from  an  independent  neutralization-reionization  ex¬ 
periment,  March  and  Young10  have  concluded  that  there  is 
insufficient  evidence  to  support  the  claim  the  H30  radical 
survives  fora  transit  time  of 0.41  /rs.  The  trace  mass  spectro- 
metric  signal  at  the  mass-to-charge  ratio  of  19  appropriate 
for  H3  0  is  identified  as  DHO  + . 

The  most  recent  theoretical  study  by  Niblaeus,  Roos, 
and  Siegbahn  (NRS)6  found  a  local  minimum  on  the  lowest 
potential  surface  of  HsO  in  geometry  20.5  keal/mol 
above  H20  +  H  with  a  barrier  to  dissociation  reported  as 
3.4  keal/mol,  which  was  expected  to  be  too  shallow  to  con¬ 
tain  a  vibrational  level.  Thus,  the  existence  of  metastable  H3 
O  in  the  gas  phase  with  a  measurable  lifetime  was  judged  to 
be  improbable. 

The  H30  molecule  is  included  in  the  recent  discussion 
by  Herzberg"  of  hypervalent  hydrides  which  might  be  ex¬ 
pected  to  exist  as  Rydberg  molecules  chiefly  because  of  the 
stability  of  the  corresponding  ground  state  ion.  However,  in 
that  review  he  reports  that  no  discrete  spectral  features  have 
been  found  and  the  broad  continuous  bands  that  were  found 
are  difficult  or  impossible  to  identify.  The  locations  of  Ryd¬ 
berg  energy  levels  of  H30  have  been  predicted  theoretically 
by  Raynor  and  Herschbach12  using  a  SCF  and  Koopmans’ 
theorem  approach. 

Finally,  the  importance  of  the  OH  +  H2  reaction  in  hy¬ 
drogen/oxygen  flames  has  led  to  a  number  of  studies13'”  of 
the  ground  state  potential  surface,  all  of  which  lead  to  a 
planar  transition  state  with  a  classical  barrier  height  of  —6 
keal/mol.  Although  at  nearly  the  same  energy,  this  transi¬ 
tion  state  is  at  a  completely  different  geometry  and  lias  a 
completely  different  character  than  the  minimum  identified 
by  Niblaeus,  Roos,  and  Sieghbahn  (NRS).4The  long  range 
part  of  the  OH-H2  potential  has  been  considered  by  Ko- 
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chanski  and  Flower14  who  were  interested  in  rotational  exci¬ 
tation  of  OH  by  Hj  at  interstellar  temperatures. 

In  this  work,  we  have  used  a  MCSCF/CI  approach  to 
survey  the  low-lying  states  of  H3  O,  paying  particular  atten¬ 
tion  to  the  symmetric  geometries  expected  for  H3+  -0“  in¬ 
teraction.  Local  minima  on  the  lowest  two  doublet  surfaces 
of  H3  O  have  been  identified  and  characterized.  The  correla¬ 
tion  diagram  linking  these  minima  to  the  different  possible 
dissociation  products  has  been  established  and  the  barriers 
to  dissociation  to  H20  +  H  and  to  OH  +  H2  have  been  de¬ 
termined  by  MCSCF  gradient  techniques.  Computational 
issues  addressed  during  the  course  of  this  study  are  also  dis¬ 
cussed. 


CALCULATIONS 

Two  Gaussian  basis  sets  have  been  used  in  our 
CASSCF/CI  calculations  of  H30.  The  initial  survey  of  the 
low-lying  potential  surfaces  of  H3C>  used  a  double  zeta  plus 
polarization 17  basis  set  augmented  by  diffuse  s  (0.08 )  and  p 
(0.06)  functions  on  oxygen,4  0(5r3pld)  H(2rlp),  denoted 
DZPR.  The  diffuse  functions  have  been  added  to  describe 
the  negative  ion  for  those  states  that  can  be  thought  of  as 
H3*  +  0“  or  to  describe  Rydberg  character.  For  the  final 
calculations  of  the  stationary  points,  the  basis  set  denoted 
DZP2R  included  an  additional  diffuse  r(0.03)  and  p(0.Q2) 
function  on  oxygen  for  more  accurate  description  of  the 
Rydberg  states. 

In  the  survey  calculations,  molecular  orbitals  optimized 
on  the  lowest  state  ( 1  *A ')  were  obtained  by  the  complete 
active  space  (CAS)  MCSCF  procedure  in  which  nine  elec¬ 
trons  were  correlated  in  a  seven  orbital  active  space.  Only  the 
orbital  arising  from  the  oxygen  lx  was  kept  frozen.  Energies 
were  obtained  from  first-order  Cl  (FOCI)  calculations  that 
included  all  single  excitations  from  the  active  space.  Except 
as  noted,  calculations  were  carried  out  in  C,  symmetry,  in 
which  the  symmetry  plane  is  perpendicular  to  the  plane  of 
the  three  H  atoms  and  includes  the  O  atom.  In  lA  '  symmetry, 
for  example,  the  MCSCF  expansion  (5a’,  2a*  active  space) 
results  in  250  configurations  and  the  FOCI  in  14  902 
( 17  080)  configurations  with  the  DZPR  and  DZP2R  basis 
sets. 

Geometries  of  the  minima  and  transition  states  were  op¬ 
timized  at  the  MCSCF  level  (DZP2R  basis)  using  analytic 
second  derivatives.  At  these  points,  energies  were  obtained 
by  multireference  single  and  double  excitation  Cl 
(MRSDCI)  calculations  using  MCSCF  orbitals  optimized 
on  the  lowest  state  of  each  molecular  symmetry.  Configura¬ 
tions  which  differed  by  at  most  two  electrons  from  a  set  of  2- 
9  reference  configurations  were  included  in  the  Cl  expan¬ 
sions,  which  totaled  —25  000-118  000  configurations.  The 
weight  of  the  reference  configurations  in  the  resulting  Cl 
wave  function  was  -0.94-0.95  for  all  calculations. 

As  discussed  below,  at  certain  geometries,  the  second 
electronic  state  is  a  3p  Rydberg  state,  i.e.,  an  H30+  core 
surrounded  by  an  oxygen  3p  Rydberg  electron.  For  Cl  calcu¬ 
lations  on  this  state  in  iA  '  symmetry,  an  0(  3p)  orbital  deter¬ 
mined  as  the  virtual  orbital  of  a  SCF  calculation  on  H,0  + 
was  added  to  the  active  space,  resulting  in  a  (6,2)  active 
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space.  The  most  similar  orbital  in  the  external  space  was 
deleted. 

A  limited  number  of  calculations  at  the  complete  sec¬ 
ond-order  Cl  (SOCI)  level  which  included  all  single  and 
double  excitationr  from  the  active  to  the  external  space  with¬ 
out  iimitation  by  reference  configurations  totalling 
—330000  configurations  with  the  DZPR  basis  were  also 
performed.  These  calculations  and  other  technical  problems 
addressed  in  the  course  of  our  study  of  H3  O  are  reported  in  a 
separate  section  below. 

Some  preliminary  geometry  optimizations  were  carried 
out  with  the  GAMESS  program1*;  final  geometry  determi¬ 
nations19  and  all  other  calculations  used  the  alchemy  it 
program  system.20 

POTENTIAL  SURFACES  AND  MINIMA 
Survey  of  Cir  geometries 

We  first  present  a  survey  of  potential  surfaces  of  H,0 
calculated  with  the  DZPR  basis  set  and  restricted  to  C3, 
symmetry.  This  restriction  is  consistent  with  our  model  of 
the  ion  pair  as  a  H3+  ion  and  is  also  consistent  with  the  sym¬ 
metric  pyramidal  geometry  optimized  by  NRS4  for  the  low¬ 
est  state.  The  geometries  have  been  characterized  in  terms  of 
two  parameters:  a,  the  H-H  distance,  i.e.,  the  side  of  the 
hydrogen  equilateral  triangle,  and  d,  the  vertical  distance 
from  the  oxygen  atom  to  the  center  of  the  triangle.  The  low¬ 
est  state  belongs  to  the  totally  symmetric  L4,  symmetry  in 

(3A '  in  C,)  symmetry.  The  second  state  is  a  doubly  de¬ 
generate  3E state  in  CjD  which  has  two  components,  the  2 2 A 
and  1 3  A  *  states  in  C,  symmetry. 

Results  for  the  lowest  potential  energy  surface  of  H  3  O  at 
various  fixed  H-H  distances  between  1.65  and  3.50  a„  are 
shown  in  Fig.  1  as  a  function  of  the  vertical  oxygen  distance. 
At  H-H  distances  near  1.65  a0 ,  the  equilibrium  separation  of 
H3* ,  the  wave  function  is  a  H3*  -O-  ion  pair  as  anticipated. 
However  the  lowest  energy  is  found  for  a  much  larger  H-h 
distance  and  closer  approach  by  oxygen,  i.e.,  a  =  3.05  a0  and 
d  =  0.92  a0,  where  the  system  can  be  described  as  a  H3  O  * 
ion  with  a  Rydberg  s  electron  on  oxygen.  This  is  in  qualita¬ 
tive  agreement  with  the  results  of  the  previous  study  by 
NRS. 

Cuts  through  the  potential  surface  at  constant  H-H  dis¬ 
tance  for  the  I2£  excited  state  are  shown  in  Fig.  2.  The  re¬ 
quired  equivalence  between  the  2  3 A  '  and  the  1 3A  *  states 
was  confirmed  in  our  calculations  in  C,  symmetry.  There  are 
two  distinct  minima  in  the  1 3E  state  in  C3l.  geometries  and 
there  is  no  p3th  without  an  energy  barrier  for  conversion  of 
the  upper  local  minimum  to  the  lower  one  in  symmetric  geo¬ 
metries.  The  upper  minimum  occurs  at  a  H-H  separation  of 
1.65  Oq.  At  this  geometry,  which  corresponds  to  the  equilib¬ 
rium  geometry  of  H3+ ,  on  the  basis  of  population  analysis  of 
the  natural  orbitals,  both  the  first  and  second  2 A  states  can 
be  identified  a5  true  ion  pairs.  (Because  of  the  open  shell 
nature  of  O-,  lsJ2sJ2 ps,  combinations  of  H3*  and  0“  can 
lead  to  more  than  one  state.)  In  these  calculations,  with  the 
DZPR  basis  set,  at  the  geometry  of  the  lower  minimum,  the 
2  3 A  'state  is  a  valence  state.  However,  as  reported  below, 
^when  provision  is  made  to  treat  the  3p  Rydberg  state  proper- 
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FIG.  I.  Cuts  through  the!  in  Cw  symme¬ 
try.  ( I  '  in  C, )  potential  surface  at  fixed  H- 
H  distance  a.  as  a  function  of  the  distance  of 
the  O  from  the  H,  plane  d,  calculated  at  the 
FOCI  level  ( DZPR  basis  set). 


d(3o) 


Iy,  our  final  calculations  indicate  the  3 p  Rydberg  state  is 
lower  in  energy  than  this  valence  minimum. 

The  potential  energy  surfaces  for  the  corresponding 
states  of  quartet  symmetry  have  also  been  investigated  and 
are  illustrated  in  Fig.  3.  As  for  the  doublet  states,  the  lowest 
state  is  the  1  *At  and  the  second  state  is  the  1  ^.correspond¬ 
ing  to  the  2  *A  and  1  *A  *  states  in  C,  symmetry.  There  is  no 
binding  in  states  of  quartet  multiplicity  in  C»  geometries. 
These  states  have  not  been  considered  further  in  this  study. 


Optimized  geometries 

From  this  survey,  there  are  three  local  minima  in  C,„ 
restricted  geometries  to  be  characterized:  the  lowest  point  on 
the  1  Vl,  (1 2A  ')  surface  and  the  two  local  minima  on  the 
1 1E  (2 2A',  1 2A  *)  surface.  When  the  C}„  symmetry  is  re¬ 
laxed,  of  course,  the  1  2A  *  and  2  2A '  surfaces  are  no  longer 
degerirate.  Optimized  geometries  in  lower  symmetry  are 
summarized  in  Table  I  and  illustrated  in  Fig.  4. 


0  2  4  6  8  10 


FIG.  2.  Cuts  through  the  I  ’£in  C,.  symme¬ 
try,  (2!A\  1  lA'  in  C,)  potential  surface  at 
fixed  H-H  distance  a,  as  a  function  of  the  dis¬ 
tance  of  the  O  from  the  H,  plane  d,  calculated 
at  the  FOCI  level  (DZPR  basis  set ). 
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FIG.  3.  Cuts  through  the  (a)  I\t,.  (b)  1  '£ 
potential  surfaces  at  fixtd  H-H  distance  a.  as 
a  function  of  the  distance  of  the  O  from  the  H, 
plane  d,  calculated  at  the  FOCI  level  ( DZPR 
basis  set).  Dashed  line  indicates  gap  between 
calculations  started  from  small  and  large  dis¬ 
tanced.  (See  the  text). 


As  determined  by  MCSCF  analytic  second  derivative 
optimization  with  the  DZP2R  basis  set,  the  true  minimum 
(all  real  frequencies)  on  the  1  2A  *  surface  has  Cj0  symmetry. 
The  MCSCF  optimized  geometry  of  the  l 3 A  '  state  obtained 
in  this  work  is  in  good  agreement  with  the  unrestricted  Har- 
tree-Fock  (UHF)  optimized  geometry  of  NRS  as  expressed 
in  Table  II  in  terms  of  the  OH  bond  length,  1.012  A  in  the 
present  work  vs  0.989  A,  NRS,  and  the  HOH  bond  angle. 
The  slightly  smaller  bond  angle,  105.1,  present  work,  vs 
109.5,  NRS,  corresponds  to  slightly  larger  vertical  distance 
(0.404  vs  0,330  A)  between  the  oxygen  and  the  plane  of  the 
three  hydrogens.  The  HH  distances  are  extremely  close  in 
the  two  optimizations.  At  the  MCSCF  level,  the  difference  in 
energy  between  the  present  structure  and  that  of  NRS  is 
1  kcal/mol.  Because  the  1 2A  ’  state  is  a  Rydberg  state,  it  is 
expected  to  have  a  very  similar  geometry  to  the  Hj  O  * 
ion,31,33  which  is  also  demonstrated  in  Table  II. 

MCSCF  gradient  optimization  predicts  two  minima 
(with  six  re;  ’  frequencies)  on  the  1  2 A  '  potential  surface. 
The  lowest  i.  .nergy  is  a  3 p  Rydberg  state  and  is  very  similar 
m  geometry  toihel  'A  'state,  the  3s  Rydberg  state.  Optimiz¬ 
ation  of  the  *A  *  state  in  C,  symmetry  will  inevitably  break 
the  CJv  symmetry.  The  geometry  difference  is  not  signifi¬ 
cant,  the  illustrations  coui  J  not  be  distinguished  in  Fig.  4.  At 


the  higher  minimum,  which  is  quite  distorted  from  C}1  sym 
metry,  the  1  2  A  '  state  is  a  valence  state.  Finally,  the  opti¬ 
mized  geometry  of  the  ion-pair  minimum  on  the  1  2E  poten 
tial  surface,  the  upper  minimum  in  Fig.  2,  has  been 
determined  within  C3„  symmetry  by  variation  of  the  param 
eters  a  and  d.  However,  the  ion-pair  state  is  not  stable  w.h 
respect  to  distortion  from  C3„  symmetry  as  reported  below 

Energies 

Total  energies  for  the  optimized  minima  from  MRSDCI 
calculations  (DZP2R  basis  set)  are  also  listed  in  Table  I 
Relative  energies  are  given  with  respect  to  the  H20  +  H 
asymptote.  The  1 3 A  '  state,  the  3s  Rydberg  state,  is  predicted 
to  lie  17.93  kcal/ir.ol  above  H20  +  H,  in  reasonable  agree 
ment  with  the  20.5  kcal/mol  value  obtained  previously  6  At 
this  geometry  the  2A '  symmetry  MRSDCI  calculation  was 
carried  out  with  the  larger  active  space  to  treat  the  3 p  Ryd¬ 
berg  state.  From  analysis  of  the  natural  orbitals  and  the  Cl 
coefficients,  the  2 3n  '  state  may  be  characterized  as  a  3 p  Ryd 
berg  state  with,  however,  some  admixture  of  valence  charac 
ter.  In  calculations  on  the  isoelectronic  system,  FH2,  Pctsa- 
lakis  et  a/.3J  found  the  first  excited  state,  the  1 2B2 ,  to  have  a 
mixed  Rydberg  and  valence  antibonding  character  The  re! 
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TABLE  I.  Optimized  geometries  ( a0 ).  total  energies,  zero-point  energies,  and  relative  energies  for  minima  on  first  two  doublet  surfaces  of  H,0. 


Total  energy 
(hartrees) 

Relative 

Zero-pt 

Energy 

(kcal/mol) 

Relative 
energy* 
with  zero  pt 
(kcal/mol) 

State 

Symmetry 

Geometry 

MCSCF  MRSDCI 

(kcal/mol) 

1  :A,  O  '/O 

C,(C,) 

O-H  t.912 

H-H  3.036 

d  0.764 

76.556  186  -76.727  647 

1793 

18  08 

23.14 

t  :£(2  *A’,t  'A')  C,,(C, ) 

O-H  1.912 

H-H  3.036 

d  0.764 

-  76.641  597 

71.94 

l  vC 

(3p  Rydberg) 

c, 

O-Hl  1.864 

0-H2  1.877 

HI-H2  3.007 

H2-H3  3.110 

76.479  853  -  76.638  687 

73.77 

23.49 

84.39 

1  :-r 

valence 

c, 

O-Hl  1.850 

0-H2  2.462 

HI-H2  3.291 

H2-H3  3.465 

76.465  506  -  76.627  716 

80  65 

17  36 

85  14 

1  "’£(2  2 A '.  1  'A 
(ion  pair) 

')C„ 

O-H  3  627 

H-H  1.650 

d  3.5 

■  ■'6.349011  -  76.531298* 

141.57* 

'  With  respect  to  H-O  +  H  MCSCF 

-  76.603  067  MRSDCI 

-76.756  204. 

*SOCI/DZPR  HjO  +  H  asymptote  =  -  76.736  837. 


ati  ve  energy  of  the  3 p  state  of  H 3  O  is  7 1 .94  or  54.0  kcal/mol 
above  the  3s  Rydberg  state.  Note  that  essentially  the  same 
prediction,  a  relative  energy  of  73.77  kcal/mol  is  given  for 
the  3 p  Rydberg  state  by  MRSDCI  calculation  on  the  2A  " 
symmetry  with  the  (5,2)  active  space  at  the  geometry  ob¬ 
tained  by  lA  '  MCSCF  optimization. 

The  non-Rydberg  minimum  on  the  1  2A  ’  potential  sur¬ 
face  is  somewhat  higher  at  80.65  kcal/mol.  Finally,  the  ion- 
pair  local  minimum  on  the  1 2E  surface  is  much  higher  in 
energy,  141.57  kcal/mol.  This  energy  was  obtained  from  a 
second-order  Cl  calculation  with  the  DZPR  basis  set, 
which,  as  shown  below,  should  be  entirely  consistent  with 
the  other  energies. 

CORRELATION  DIAGRAM 

The  correlation  diagram  connecting  the  1  }At  and  l  !E 
states  of  HjO  and  the  corresponding  states  in  C,  symmetry 
to  the  states  of  the  possible  dissociation  limits,  OH  4-  H:, 
HjO  +  H,  Hj  +  O,  K/  4-  Q~»  and  0  +  H  +  H2»  based 
only  on  symmetry  considerations  is  shown  in  Fig.  5.  Ener¬ 
gies  of  the  asymptotes  with  respect  to  H20  +  H,  excluding 
zero  point  energy,  summarized  in  Table  HI,  have  been  ob¬ 
tained  as  follows:  The  experimental  heats  of  formation,” 
modified  by  the  experimental  zero  point  energies”-”  are  tak¬ 
en  for  the  OH(*  Jfl)  +  H2  and  the  H2  4-  H  +  O  asymp¬ 
totes.  Spectroscopic  T,  values  have  been  used  for  the  excited 
states  of  OH  and  H20.  The  ground  state  surface  of  Ht  is 
repulsive  with  respect  to  H2  4-  H.  The  H3  +  O  asymptote 


has  been  placed  according  to  the  second  order  Cl  calculated 
energy  difference  between  H  at  the  H,"  equilibrium  geome¬ 
try,  a  —  1.65  a0,  and  that  of  H2  4-  H.  Finally,  the  accurate 
calculated  value  for  the  Hj*  dissociation  energy  7  and  the 
experimental  oxygen  electron  affinity5  have  been  taken  for 
the  ion-pair  limit.  The  correlations  drawn  for  the  1  2E  state 
on  the  basis  of  symmetry  apply  as  well  to  the  ion-pair  mini¬ 
mum  (not  shown)  which  lies  ~  70  kcal/mol  higher  in  ener¬ 
gy- 

The  H,  4-  Oand  Hs*  4-  O"  limits,  which  are  very  high 
in  energy,  each  give  rise  to  an  A ,  and  E state  and  thus  corre¬ 
late  in  symmetry  to  the  lowest  states  of  H,0.  Of  greater 
interest  are  the  correlations  to  the  H20  +  H  and  OH  +  H2 
asymptotes  which  have  been  drawn  in  C,  symmetry.  Using 
the  C,  state  designations,  the  1  2A '  state  correlates  with  the 
ground  state  of  OH  +  H2  and  of  H20  +  H  and  is  higher  in 
energy  than  either  of  these  limits.  The  2 2 A '  state  is  bound 
with  respect  to  the  excited  states  of  OH  and  H20,  with 
which  it  correlates.  The  1 2A  "  state,  however,  correlates  to 
the  OH  ground  state,  which  has  both  an  A  ‘  and  A  '  compo¬ 
nent.  In  the  other  direction,  the  !  2 A  “  state  correlates  to  a 
highly  excited  state  of  H2Q  +  H.  Coupling  between  the 
1 2A  '  and  2 2 A '  states  would  cause  the  stability  of  the  2 2 A ' 
state  to  depend  on  that  of  the  1 2A  '  state.  Anticipating  our 
results  for  the  ion-pair  local  minimum,  we  also  note  that  the 
!  2E,  state  correlates  to  the  H2  4-  H  +  0  limit  in  C,  (no) 
symmetry. 

From  this  correlation  diagram,  barriers  in  the  dissocia- 
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FIG.  3.  Correlation  diagram  linking  the  I  *A,  and  the  I  Estates  of  H,0  to 
asymptotic  limits  on  the  basis  of  symmetry.  The  H,*  +  O"  and  H,  +  O 
limits  have  been  connected  in  C,.  symmetry,  the  OH  +  H:  and  H2  O  +  H 
limits  have  been  connected  in  C,  symmetry  and  the  H,  +  H  +  O  limit  in  C, 
(no)  symmetry. 


FIG  4.  Geometries  optimized  by  MCSCF  gradient  calculations  (DZP2R 
basis)  (a)  1  ‘A 'minimum. (b)  1 XA  * 3pRydbergrmntmum,(c)  1  3  A  'mini¬ 
mum  of  valence  character,  (d)  transition  state  1  ‘A  -  H20  ■+•  H,  (e)  Sta- 
uonary  point  with  two  imaginary  frequencies  determined  in  search  for  tran¬ 
sition  state  1  ‘A  '-OH  +  H,  (see  the  text),  and  (0  transition  state  1  !A  * 
valence  minimum— OH  +  H  +■  H. 

tion  of  the  1 2A '  state  to  H20  +  H  and  to  OH  +  H2  and  of 
the  1 2 A  '  state  to  OH  +  H2  need  to  be  determined.  For  the 
1  lA  *  state,  the  relationship  between  the  3 p  Rydberg  mini¬ 
mum  and  the  valence  minimum  also  needs  to  be  investigat¬ 
ed. 

STABILITY  OF  POTENTIAL  MINIMA:  DISSOCIATION 
BARRIERS  AND  RADIATIVE  LIFETIMES 

Dissociation  of  the  1  * A  state 

The  transition  states  for  dissociation  of  the  minimum  on 
the  lowest  potential  energy  surface  to  H20  +  H  and  to 


TABLE  II.  Comparison  of  H,0  I  * A  ’  and  HjO*  geometry  (bond  length  in 
A,  angles  in  degrees)  and  zero-point  energy  (kcal/mol). 


H,0(!  3A') 

H,0* 

Present  work  Ref.  6* 

Calculated* 

Exptr 

O-H 

1.012  0.989 

0.987 

zHOH 

Zero- point 

103.1  109.3 

1 1 1.6 

111.3 

energy 

18.08 

20.32 

*UHF  optimization. 

‘DeFreesand  McLean,  Ref.  21.  6-3 IG*  basis  Cl  opt. 

'Symons.  Ref.  22  deduced  from  nuclear  magnetic  resonance  (NMR)  data. 


OH  -f  Hj  have  been  investigated  by  MCSCF  gradient  opti¬ 
mization.  Geometries,  energies,  and  zero  point  energies  ( 5 
or  4  real  frequencies)  of  the  stationary  points  are  listed  in 
Table  IV.  A  very  low  barrier,  3.58  kcal/mol,  is  found  for 
dissociation  to  HzO  -f  H,  which  has  a  conventional  transi¬ 
tion  state  with  one  imaginary  frequency  (23 10/).  The  geom¬ 
etry  of  the  transition  state  is  in  between  that  of  the  1  2 A  ' 
minimum  and  of  H20.  The  symmetric  hydrogens,  H2  and 
K3  in  Fig.  4,  form  the  H20  molecule.  Including  the  zero- 
point  energy  which  is  larger  by  3.2  kcal/mol  for  the  1  2A 
equilibrium  than  for  the  transition  state,  the  barrier  height  is 
reduced  to  0.40  kcal/mol.  This  value  is  smaller  than  uncer¬ 
tainties  in  the  calculation.  Note  that  the  transition  state  ge- 


TABLE  III.  Energies  of  asymptotes  in  kcal/mol  relative  to  H.O  +  H  ex¬ 
cluding  zero- point  energy. 


Asymptote 

Relative  energy 

HjOUU.J+H 

0.0 

HjO  (A  ’B,)  +  H 

133.71* 

oh  (jr’m  +  Hj 

160* 

OH  (A  2i*)  +  H2 

10948' 

H,  +  H  +  O 

122.67* 

H,  +  0 

204  57* 

H,*  +0- 

296.30* 

'Herzberg  (Ref.  26). 

‘Heats  of  formation  ( Ref  24)  modified  by  zero  point  energies  ( Refs.  23  and 
26). 

'Huber  and  Herzberg  (Ref.  23). 

*0x10018164  (present  work). 

'H/  dissociation  energy  Dykstra  el  a!.  (Ref.  27).  O  electron  affinity  !  Ref. 
5). 
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TABLE  IV  Optimized  geometries  (o0).  total  energies,  zero-point  energies,  and  relative  energies  for  transition  states  and  stationary  points. 


Total  energy 
(hart  reel) 

Relative 

Zero-pt 

Transition 

Geometry 

MCSCF 

MRSDCI 

energy 

(keal/mo!) 

(keal/mol) 

(keal/mol) 

w/zero  pt 

V4-H.O  +  H 

O-HI  2.209* 
0-H2  1.877 

H1-H2  3.245 
H2-H3  2.969 

-76.552494 

-76.721  953 

21.50 

14.90 

3  58 

0.40 

1  ^'—OH  +  H  +  H 

O-Hl  I.8JS 
C-H2  4.241 
H1-H2  4.675 
H2-H3  4.014 

-76.434  859 

-  76.572  264 

115.46 

10.65 

34.81 

28.09 

1 2 A  '  Rydberg  - 
1  2 A  '  valence 

O-Hl  1.844 
0-H2  2.198 
H1-H2  3.153 
H2-H3  3.264 

-76.463  780 

-76  632432 

77  69 

16.79 

Stationary  point* 

O-HI  1.865 
0-H2  2.373 
H1-H2  3.247 
H2-H3  2.131 

-76.515  968 

-  76.687  349 

43.22 

i 

*t 

. 

i 

11  38 

25.30 

18  60 

“With  respect  to  HjO  +  H  MCSCF  -  76.603  067  MRSDCI  -  76.736  204.  I 

6NRS  (Ref.  6)  geometry  O-Hl  =  2.350, 0-H2  =  1.821,  H1-H2  =  3.422,  H2-H3  =  2.943. 

‘Stationary  point  with  two  imaginary  frequencies  located  in  search  for  transition  state  1  2 A  '-OH  +  H2  (see  the  text). 


ometry  is  quite  similar  to  that  found  by  NRS  who  reported  a  l 2 A  ’  state  depends,  in  any  event,  on  the  very  low  barrier  to 

barrier  height  of  3.4  keal/mol  without  zero-point  energy,  in  formation  of  H2  O  +  H. 
excellent  agreement  with  present  results.  We  also  find,  in 

agreement  with  NRS,  that  the  change  from  Rydberg  to  va-  Dissociation  of  th®  1 7 A  '  state 

lence  character  as  the  1 2A '  state  dissociates  is  due  to  a  The  dissociation  path  of  the  lowest  potential  surface  of 

change  in  the  molecular  orbital  and  not  to  configuration  2A  '  symmetry  is  more  complicated.  As  shown  in  Table  IV, 

mixing,  an  example  of  Mulliken’s  **MO  Rydbergization."2*  two  transition  states  have  been  identified  at  the  MCSCF  lev- 

A  similar  situation  has  been  discussed  in  the  dissociation  of  el.  The  one  listed  at  ;he  bottom  of  Table  IV  connects  the 

NH«  to  NHj  +  H.w  Rydberg  and  valence  minima  on  the  1  2A  '  surface.  The  ge- 

We  were  not  successful  in  locating  the  transition  state  ometry  is  in  between  that  of  the  two  minima.  It  has  an  imagi- 

for  dissociation  of  the  1  2 A '  state  to  OH  +  H2 .  The  station-  nary  frequency  of  1 672 /  and  distortion  of  the  transition  state 

ary  point  characterized  by  two  imaginary  frequencies,  3493/  geometry  in  the  two  directions  according  to  the  imaginary 
and  1520/,  that  has  been  determined  is  reported  in  Table  IV.  mode  leads  to  the  two  minima.  The  transition  state  lies  ap- 
In  our  search  for  the  transition  state,  the  stationary  point  proximately  1  keal/mol  above  the  1 3 A  “  valence  minimum  at 

geometry  was  distorted  by  varying  amounts  according  to  the  the  MCSCF  level.  However,  MRSDCI  calculations  place 

normal  modes  associated  with  each  imaginary  frequency,  in  the  transition  state  2.96  keal/mol  below  the  energy  of  the 

turn,  to  provide  a  starting  point  for  further  optimization.  We  valence  minimum.  Thus,  the  true  status  of  the  valence  mini- 
found,  however,  that  distortion  of  the  transition  state  geom-  mum  as  a  stationary  point  on  the  1  2 A  '  potential  surface  is  in 

etry  in  one  direction  according  to  the  normal  mode  corre-  doubt. 

sponding  to  the  smaller  frequency,  which  preserves  C,  sym-  The  other  transition  state  on  the 2 A  *  surface  lies  at  rath- 

metry,  led  to  optimization  to  OH  -f  H2 ,  and  distortion  in  the  er  high  energy,  ~35  keal/mol  above  the  valence  minimum, 

other  direction  led  back  to  the  1  2A  minimum.  Distortion  of  The  transition  state  has  a  single  imaginary  frequency,  1854/. 

the  transition  state  according  to  the  normal  mode  of  the  larg-  However,  at  the  transition  state  geometry,  the  H2  internu- 

er  imaginary  frequency,  which  breaks  C,  symmetry,  led  to  clear  separation  is  extremely  large,  4.0  a0,  and  the  MCSCF 

formation  of  HjO  +  H  with  one  of  the  symmetric  hydrogens  energy  is  within  1  keal/moi  of  that  ofOH  +  H  +  HwithOH 

becoming  the  hydrogen  atom  asymptotically.  Distortion  in  at  its  equilibrium  intemuclear  separation.  Thus,  it  appears 

the  opposite  direction  led  to  elimination  of  the  other  sym-  that  the  transition  state  that  has  been  identified  leads  to 

metric  hydrogen.  The  energy  ofthe  stationary  point  is  25.30  OH  +  H  +  H  and  not  to  OH(A,2ID  +  H2  as  anticipated 

keal/mol  above  that  of  the  1 2A  minimum.  Stability  of  the  from  the  correlation  diagram  in  Fig.  5.  Distortion  according 
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to  the  imaginary  mode  in  the  other  direction  definitely  leads 
back  to  the  1 2 A  ’  valence  minimum.  It  was  not  possible  to 
locate  a  transition  state  with  smaller  H-H  distance  within  C, 
symmetry. 

From  these  results,  it  is  not  possible  to  unambiguously 
describe  the  dissociation  of  the  l2A  *  state.  Further  explora¬ 
tion  at  the  Cl  level  would  be  required  to  determine  if  a  va¬ 
lence  minimum  exists.  More  imporantly,  investigation  in 
lower  symmetry  would  be  required  to  search  for  a  lower 
energy  path  to  OH(-SfJn)  +  H2 

Dissociation  of  tht  ion-pair  minimum 

The  ion-pair  minimum  on  the  1  2E  potential  energy  sur¬ 
face  lies  higher  in  energy  than  the  H2  +  H  +  O  asymptote 
and  is  connected  to  it  by  symmetry.  The  plot  in  Fig.  6  dem¬ 
onstrates  that  distortion  of  the  Cu  geometry  of  this  local 
minimum  leads  to  dissociation  of  the  2  2 A  ’  state  (22A  in  C, 
symmetry)  to  H2  +  H  -H  O  without  a  barrier.  The  dissocia¬ 
tion  path  chosen  is  as  follows.  FOCI  energies  obtained  with 
the  DZPR  basis  set  have  been  plotted.  Starting  from  the 
equilibrium  geometry,  a  —  1.65  a0,  d  =  3.5  a0  and  keeping 
the  vertical  distance,  d,  and  the  two  sides  of  the  triangle 
constant,  the  triangle  has  been  opened  gradually  with  a  base 
going  from  1.65  to  3.30  a0  where  the  three  H’s  are  aligned 
and  kept  equidistant  from  the  central  H.  Then  one  of  the  H’s 
and  the  O  were  gradually  pulled  out  while  the  two  remaining 
hydrogens  were  moved  closer  until  the  equilibrium  geome¬ 
try  of  H2  (1.402  a0)  was  attained.  The  reaction  coordinate 
reported  in  Fig.  6  is  the  distance  between  the  two  most  dis¬ 
tant  hydrogens.  While  this  may  not  necessarily  provide  the 
minimum  energy  path,  it  is  sufficient  to  establish  that  the 
ion-pair  state  of  H,  O  is  not  expected  to  be  stable. 

The  3  2A  state  ( 1  2A  '  in  C,  symmetry),  also  plotted  in 


Fig.  6,  shows  a  barrier  on  the  chosen  dissociation  path  due  to 
an  interaction  with  a  higher  ion-pair  state.  We  have  not  opti¬ 
mized  the  path  for  this  state  because,  in  any  event,  coupling 
with  the 2 A '  state  should  prevent  stability  of  the 2  A  ‘  compo¬ 
nent  of  the  ion-pair  minimum. 

Radiative  lifetimes 

Transition  dipole  matrix  elements  connecting  the  2  2 A  ' 
and  the  l2A°  states  with  the  1 14'state  have  been  calculated 
from  FOCI  wave  functions  at  the  geometries  of  the  upper 
state  minima  and  used  to  estimate  oscillator  strengths  and 
radiative  lifetimes,  listed  in  Table  V,  from  the  very  approxi¬ 
mate  formulas  given  there.  The  3 p  Rydberg  state  has  an  esti¬ 
mated  lifetime  of  18  ns  which  corresponds  to  a  transition 
probability  of 5.46  X 107  s"  ’.  Raynor  et  al.n  obtained  a  tran¬ 
sition  probability  of 5.64x  107  for  this  state.  Note  that  essen¬ 
tially  the  same  result  is  obtained  for  the  3 p  state  whether 
calculated  in  2 A '  or  2A  '  symmetry.  At  the  geometry  of  the 
l2 A  '  valence  minimum,  the  calculated  transition  moment  is 
approximately  half  as  large  which  would  correspond  to  a 
radiative  lifetime  of  34  ns  if  the  upper  state  is  considered  as  a 
bound  state. 

COMPUTATIONAL  ISSUES 

Although  the  final  energies  at  the  stationary  points  were 
determined  by  MRSDCI  calculations,  a  great  deal  of  explo¬ 
ration  was  performed  at  the  FOCI  level.  The  comparisons  of 
total  and  relative  energies  at  the  1  2 A  '  equilibrium  geometry 
and  at  the  H2  O  +  H  and  OH  +  H2  asymptotes  at  their  spec¬ 
troscopic  geometries  reported  in  Table  VI  provide  the  cali¬ 
bration  of  FOCI  results. 

The  most  important  comparison  is  that  between  the 


FIG.  6.  Dissociation  of  ihc  upper  local  mini¬ 
mum  on  the  1  ‘E  potential  surface  of  H ,  O  to 
Hj  +  H  +  O  as  a  function  of  the  pseudo  reac¬ 
tion  coordinate  5  defined  in  the  figure,  calcu¬ 
lated  at  the  FOCI  level  (DZPR  basis  set). 
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TABLE  V.  Transition  dipole  matrix  elements,  Mediator  strengths  and  esti¬ 
mated  lifetimes  for  2  lA '  and  1  lA  *  states  of  H  30. 


Transition  Geometry 

fiTT'.V 

/1(A)* 

r 

r(ns)‘ 

21A'—  l*A‘  !J,4'min 

1311 

5824 

0.278 

IS 

1V4'-|H'  l'A‘3p  Rydberg 

1413 

5803 

0.303 

17 

1  1  U  '  valence 

t.112 

4382 

0.086 

34 

*  Ca’  ‘laud  from  electronic  energy  difference. 

‘E:  jatedas/^jAECS^r,)1). 

' Given  by  r-'  (*-')  -  [f/X'{k)\  6.670 246 X  10,J. 


(6,2)FOCI  and  the  (6,2)MRSDCI  at  the  1 2A '  minimum, 
both  of  which  have  a  larger  active  space  to  accomodate  the 
3 p  Rydberg  orbital.  The  relative  energy  of  the  1 2A '  state 
with  respect  to  HjO  -f  H  changes  by  Jess  than  0.2  keal/mol 
on  inclusion  of  double  excitations  in  the  Cl  expansion,  while 
the  difference  in  relative  energy  of  the  2  2 A’  state,  the  3 p  Ryd¬ 
berg  state,  is  4  kcal/moL  The  Cl  coefficients  and  natural 
orbitals  were  extremely  similar  in  the  two  calculations. 
Thus,  the  FOCI  provides  meaningful  results  for  the  H3  O 
molecule.  Comparing  the  (5.2)  and  (6,2)  FOCI  results,  we 
observe  that,  as  expected,  the  total  energy  of  the  1  2 A '  state  is 
virtually  unaffected  by  the  inclusion  of  the  3 p  orbital. 

The  transition  dipole  moment  connecting  the  1 2A  ’  and 
2 2 A '  states  is  unaffected  by  the  double  excitations.  A  value  of 
2.3  U  ea0  is  obtained  from  FOCI  calculation  as  compared 
with  2.315  ea0  from  the  MRSDCI.  Therefore,  transition  di¬ 
pole  moments  from  FOCI  calculations  with  a  common  orbi¬ 


tal  set  (2A '  MCSCF  orbitals)  were  reported  in  Table  V.  For 
calculation  of  the  3 p  Rydberg  state  in  2 A  *  symmetry,  an  a" 
symmetry  3 p  SCF  virtual  orbital  was  added  to  the  active 
space,  resulting  in  a  ( 5a',  3 a" )  active  space. 

A  few  additional  comparisons  are  also  noted  in  Table 
VI.  Limiting  the  FOCI  expansion  to  at  most  double  excita¬ 
tions  with  respect  to  the  same  set  of  reference  configurations, 
labeled  FOCI  (refs),  is  seen  to  have  a  negligible  effect,  as 
compared  with  the  full  FOCI.  This  limitation  was  imposed 
in  the  MRSDCI  calculation.  At  the  asymptotes,  the 
MRSDCI  result  is  compared  with  a  complete  second  order 
Cl  calculation  (SOCI),  which  includes  all  single  and  double 
excitations  from  the  active  space  to  the  interna]  space  with¬ 
out  restriction  by  reference  configurations.  Introduction  of 
reference  configurations  modifies  the  total  energy  by  0.6 
mhartrees  and  the  relative  energy  of  the  asymptotes  by 
<0.1kcal/mol.  The  effect  of  selecting  reference  configura¬ 
tions  has  recently  been  discussed  by  Taylor  and  Partridge. w 
The  SOCI  calculation  was  done  with  the  smaller  DZPR  ba¬ 
sis  set.  The  comparison,  however,  is  valid  because  at  the  as¬ 
ymptote,  use  of  the  DZP2R  basis  set  changes  the  MCSCF 
and  FOCI  energies  by  <0.2  mhartrees. 

The  OH  +  H2  excitation  energy  with  respect  to 
HjO  +  H  is  calculated  at  the  FOCI  level  to  be  16.13  kcal/ 
mol,  in  fortuitously  close  agreement  with  the  experimental 
result  of!  6.0  keal/mol,  determined  by  modifying  the  experi¬ 
mental  heat  of  formation  by  the  zero-point  energies.  Intro¬ 
duction  of  double  excitations  reduces  the  relative  energy  by 
4.7  keal/mol.  The  FOCI  result  is  consistent  with  the  results 
ofWalch  and  Dunning13  who  obtained  16.7  keal/mol  with  a 


TABLE  VI.  Comparison  of  calculation  at  1  ‘  minimum  and  at  HjO  +  H  and  OH  +  H,  asymptotes. 


l  2a- 

2 2  A  '* 

Total  energy 
(hanrees) 

Relative  energy* 
(keal/mol) 

Total  e.iergy 
(hanrees) 

Relative  energy' 
(keal/mol) 

(5.2)  FOCI 
(6  2)  FOCI 

-0.634  937 
-  0.636  162 

18.77 

18.01 

-  0.556  704 

49.87 

(6,2)  FOCI  (refs) 

-0.635  913 

18.10 

-  0.556  224 

50  03 

(6.2)  MRSDCI 

-  0.727  647 

17  93 

-  0.641  597 

54  01 

(6.2)  FOCI/{6,2) 
orbitals 

-  0.654  706 

HjO  +  H* 
Total  energy 
(harness) 

20.51 

Total  energy 
(hartrees) 

OH  +  Hj* 

Relative  energy* 
(keal/mol) 

(5.2)  FOCI 

-0.664  851 

-  0.639  148 

16.13 

(5.2)  FOCI  (reft) 

-  0.664  765 

-0.639  098 

16.11 

(5.2)  MRSDCI 

-0.736  204 

-0.738  047 

11.40 

(5.2)  SOCI/DZPR 

-0.756  837 

-  0.738  573 

11.46 

(6.2)  FOCI/(6,2) 

-  0.687  379 

orbitals 

Expt 

16.0 

*  1 2 A  ’  minimum  geometry. 

‘With  respect  to  HjO  +  H. 
e  With  respect  to  1  *A 

*  HjO  geometry:  OH  a  1.810  0,0  =  104.3. 
V(OH)  =  1.833  OoriHH)  -  1.402a,. 
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comparable  calculation  (POL-CI)  and  comparable  size  ba¬ 
sis  set.  Unpublished  results  of  Kraka  and  Dunning14  indicate 
a  larger  basis  set,  at  least  triple  zeta  +  polarization  is  needed 
to  obtain  an  excitation  energy  within  1  kcal/mol  of  experi¬ 
ment  when  double  excitations  are  included  in  the  wave  func¬ 
tion. 

All  of  our  calculations  use  MCSCF  orbitals  obtained 
with  a  { 5,2)  active  space.  Because,  the  1  2 A '  minimum  is  a  3s 
Rydberg  state,  at  this  geometry,  the  complete  active  space 
should  actually  include  an  additional  orbital  (6,2)  to  acco¬ 
modate  the  Rydberg  orbital.  At  the  asymptote  or  at  any  va¬ 
lence  geometry,  the  larger  active  space  would  contain  an 
“extra”  orbital.  This  is  a  persistent  problem  in  the 
MCSCF/CI  treatment  of  Rydberg  states.  The  relative  ener¬ 
gy  of  the  1 2 A '  3s  state  with  respect  to  the  H20  +  H  asymp¬ 
tote  has  been  determined  by  a  (6,2)  FOCI  using  MCSCF 
orbitals  determined  with  the  (6,2)  activespace.  As  shown  in 
Table  VI,  use  of  the  larger  active  space  changes  the  excita¬ 
tion  energy  by  only  2.5  kcal/mol.  A  previous  study31  of  NH2 
demonstrated  that  use  of  a  larger  active  space  produces  a 
tolerable  variation  in  the  interaction  energy  of  valence  geo¬ 
metries,  even  at  the  MCSCF  level.  We  may  also  note  the  3s 
Rydberg  orbital  from  the  larger  MCSCF  is  very  similar  to 
that  obtained  from  the  (5,2)  calculation.  The  additional  or¬ 
bital  in  the  (6,2)  MCSCF  is  the  expected  antibonding  com¬ 
bination  of  the  oxygen  p  and  the  H’s  to  complete  the  active 
space. 

The  FOCI  prediction  illustrated  in  Fig.  6  that  the  ion- 
pair  state  dissociates  without  a  barrier  to  H2  +  H  +  O  needs 
to  be  critically  examined  because  the  FOCI  wave  function 
does  net  give  the  correct  sign  for  the  electron  affinity  of  the 
isolated  oxygen  atom.  At  the  SOCI  level,  the  electron  affin¬ 
ity  of  oxygen  is  0.82  eV  with  the  DZPR  basis  and  0.84  eV 
with  the  DZP2R  basis  in  qualitative  agreement  with  the  re¬ 
ported  value,3 1.46  eV.  The  discrepancy  with  experiment  is, 
however,  small  compared  with  the  7.77  eV  coulombic  at¬ 
traction  at  the  2  2A '  ion-pair  local  minimum. 

In  Table  VII,  the  relative  energy  of  the  2 2 A  '  state  calcu¬ 
lated  at  the  FOCI  and  SOCI  level  (DZPR  basis)  is  com¬ 
pared  at  the  ion-pair  minimum  and  at  two  critical  points 
along  the  dissociation  path  in  Fig.  6.  Truncating  the  Cl  ex¬ 
pansion  causes  a  maximum  discrepancy  of  2.3  kcal/mol  in 
the  2  2A  '  state  confirming  the  qualitative  conclusion  that  the 
ion-pair  local  minimum  dissociates  without  a  barrier.  We 


TABLE  VII.  Relative  energy  (kcal/mol)  of  2  state  with  respect  to 
HjO  +  H. 


FOCI 

SOCI 

H2-H3  («„)' 

DZP2R* 

DZPR* 

DZPR* 

1.63 

143.53 

143.90 

141.57 

2.50 

138.54 

137.83 

3.50 

133.37 

133.99 

*HjO  +  H  energy:  —  76.664  851  hartrees. 
‘HjO  +  H  energy:  —  76.664  667  hartrees. 
*  HjO  +  H  energy:  -  76.756  837  hartrees. 
'Corresponds  to  coordinate  5 in  Fig.  6. 


note  that  accurate  calculation  of  the  electron  affinity  of  the 
oxygen  atom  remains  a  classic  problem53  in  quantum  chem¬ 
istry  requiring  very  extended  basis  sets  and  extensive  treat¬ 
ment  of  correlation.  Recent  studies  have  used /  functions*5 
and/  g,  and  h  functions!34  Such  extensions  are  beyond  the 
scope  of  the  present  study. 

Finally,  we  note  that  it  was  necessary  to  monitor  the 
MCSCF  wave  function  when  approaching  dissociation  lim¬ 
its  from  smaller  intemuclear  separations,  i.e.,  from  the  left  in 
Figs.  1-3  and  6.  Particular  difficulties  were  observed  in  the 
original  survey  calculations  at  large  H3-0  distances.  The 
wave  fuiction  remained  an  ion-pair  instead  of  attaining  the 
neutral  H3  +  O  asymptote,  which  is  lower  in  energy.  In 
mapping  the  ion-pair  dissociation,  Fig.  6,  the  H3  antibond¬ 
ing  orbital,  which  should  become  an  H2  antibonding  orbital 
at  the  H2  +  H  +  O  asymptote,  became  instead  an  oxygen  2 p 
correlating  orbital.  To  overcome  these  difficulties,  the  calcu¬ 
lations  were  started  at  long  range  from  localized  orbitals  and 
carried  inward  until  the  results  matched  those  started  from 
short  range  calculations.  The  small  region  not  cohered  in 
some  of  our  calculations  on  the  quartets  is  indicated  by  a 
dashed  line  in  Fig.  3. 

DISCUSSION 
Ion-pair  states 

While  one  of  the  motivations  for  this  work  was  to  inves¬ 
tigate  a  possible  energetic  ion-pair  state  of  H3G,  our  results 
have  shown  that  there  is  not  even  a  local  minimum  with  an 
ion-pair  configuration  on  the  lowest  potential  surface.  The 
local  minimum  on  the  excited  1  2E  surface  in  C3„  geometry 
corresponding  to  H34-0"  is  not  stable  in  lower  symmetry 
Analysis  of  the  wave  functions  at  points  along  the  dissocia¬ 
tion  path  in  Fig.  6  shows  that  on  distortion  of  the  symmetric 
geometry,  the  system  quickly  loses  its  ion-pair  character, 
becoming  neutral  H3  -O.  The  lowest  state  of  H 3  is  extremely 
repulsive  with  respect  to  H2  +  H.  At  the  geometry  of  the 
H34  equilibrium,  for  example,  the  H3  ground  state  lies  8 1  9 
kcal/mol  above  H2  +  H.  Hence  the  H3-0  system  readily 
dissociates  to  H2  +  H  +  O. 

One  may  speculate  that  in  any  system  where  charge 
transfer  leads  to  unstable  neutral  species,  ion-pair  states  may 
not  be  stable.  We  have  therefore  undertaken  calculations'5  '6 
on  the  Li3  H  system  as  a  counter  example  because  both  Li3* 
and  Li3  are  known  to  be  stable.  As  will  be  reported  else¬ 
where,34  because  the  lowest  state  of  Li3  H  lies  below  all  possi¬ 
ble  asymptotes,  it  is,  indeed,  predicted  to  be  a  stable  species. 
However,  it  is  not  particularly  energetic  and  the  wave  func¬ 
tion  involves  substantial  mixing  between  neutral  and  ion- 
pair  configurations. 

Relationship  to  other  work 

Our  principal  results  for  minima  and  transition  states 
are  displayed  graphically  in  Fig.  7.  We  have  also  included 
the  energy  of  the  H3  O  4  ion  placed  by  the  experimental  pro¬ 
ton  affinity  corrected  for  zero-point  energy.3 1  The  planar 
transition  state  for  the  reaction  OH  +  H2—H20  -f  H  would 
appear  at  —  22  kcal/mol  in  Fig.  7  according  to  the  reported 
5.1-62  kcal/mol13'13  barrier  height,  approximately  4  kcal/ 
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FIG  ’  Minimi,  transition  states,  and  stationary  points  obtained  in  this  work  (MRSDCI  energies)  m  relation  to  asymptotic  energies.  The  dashed  iiansmon 
state  on  the  i  *A  '  potential  surface  connects  to  the  OH  4-  H  +■  H  asymptote  (not  shown).  (See  the  text.) 


mol  above  the  1 2 A '  state  of  H30  There  does  not,  however, 
appear  to  be  any  connection  between  this  planar  geometry 
and  the  pyramidal  geometry  of  the  H3  O  Rydberg  states.  In  a 
brief  search,  we  were  ./table  to  find  any  low  energy  path 
connecting  the  two  geometries.  Thus,  we  do  not  believe  the 
planar  transition  state  is  relevant  to  the  following  discussion. 

The  present  results  may  be  discussed  in  view  of  the  prin¬ 
cipal  experimental  findings  for  H30  in  the  literature.  Herz- 
berg1 1  has  described  an  absence  of  discrete  spectral  features. 
Gellene  and  Porter7  have  reported  observation  of  metastable 
DjO  with  a  lifetime  >0.6  fi s  by  charge-transfer  and  colli¬ 
sions!  reionization  in  K  but  find  no  evidence  for  H3  0.  In  a 
subsequent  paper,  Raskit  and  Porter*  reported  charge- 
transfer  in  Na  leads  to  observation  of  D3 160  but  not  of 
D,  "O  under  identical  conditions.  The  work  of  March  and 
Young10  leads  to  the  conclusion  that  H3  O  is  not  stable. 

The  discussion  to  date  by  Porter  and  co-workers7  *  of  the 
existence  of  H30  has  been  based  on  the  possible  barrier  to 
dissociation  to  H20  +  H  on  the  lowest  potential  surface  of 
H3  O  assuming  that  the  lowest  state  is  being  formed. 

The  present  work  confirms  the  undemanding  presented 
previously,6  i.e.,  that  the  vibrationless  potential  has  a  very 
small  barrier,  3.58  keal/mol  from  the  present  calculations 
for  dissociation  of  the  1 2 A '  state  of  H20  +  H.  The  zero- 
point  energy  of  the  1 2 A '  state  is  calculated  to  be  18.08  kcal/ 
mol  while  that  of  the  five  real  frequencies  for  the  transition 
state  is  14  90  keal/mol,  At  the  transition  state,  not  only  is 
there  one  less  bound  mode,  but  the  remaining  vibrations  are 
slightly  lower  in  energy.  Including  the  vibrational  energies,  a 
barrier  height  of  0  4  keal/mol  is  predicted.  Scaling  the  vibra¬ 
tional  energies  at  the  minimum  and  the  barrier  by  the  square 
root  of  the  reduced  mass,  which  should  be  at  least  approxi¬ 
mately  valid,  gives  a  barrier  height  including  vibration  for 


D3  0  of  1.32  keal/mol.  Considering  that  the  uncertainties  in 
the  calculation  are  larger  than  the  barrier  heights  predicted 
for  H30  or  D30  and  that  dissociation  by  tunneling  should 
be  possible  for  both  species,  it  is  not  possible  to  predict 
whether  these  calculated  barriers  on  the  1 2A  '  surface  ac¬ 
count  for  the  observation  of  D3Ot  but  not  H30.  Further¬ 
more,  any  argument  based  on  the  spacing  of  vibrational  lev¬ 
els  in  D30  would  also  predict  D3 180  to  be  stable,  which  was 
not  observed.* 

The  dissociative  nature  of  the  lowest  state  would  pre¬ 
sumably  explain  the  failure  to  observe  any  discrete  spectra 
terminating  on  the  lowest  state.  Herzberg,  in  fact,  speculates 
that  predissociation  of  possible  upper  states  leads  to  the  lack 
of  discrete  lines.  We  have  examined  the  3 p  state,  one  likely 
upper  state,  in  the  present  work.  The  only  candidate  path  is 
dissociation  along  the  1 2A "  potential  surface  to  the 
OH(X)  -f  H2  asymptote  to  which  it  correlates  directly.  Un¬ 
fortunately,  working  in  C,  symmetry,  we  have  not  been  able 
to  locate  a  transition  state  that  leads  directly  to  this  asymp¬ 
tote.  However,  it  may  be  argued  that  because  the 
OH(i!Q  +  H2  asymptote  involves  no  Rydberg  character, 
any  path  from  the  3 p  Rydberg  minimum  to  that  asymptote 
would  have  to  go  through  a  change  to  valence  character. 
Assuming  the  valence  minimum  located  by  MCSCF  gradi¬ 
ents  is  at  least  approximately  correct,  it  provides  an  estimate 
of  the  minimum  energy  of  a  valence  state  of 2 A  "  symmetry 
and  thus  a  lower  limit  to  the  dissociation  barrier.  When  zero 
point  energy  is  included,  as  shown  m  Table  I,  the  1 2 A  '  Ryd¬ 
berg  and  valence  minima  are  at  virtually  the  same  energy  so 
that  the  3 p  Rydberg  state  may  dissociate  without  a  bamcr. 

While  our  results  for  dissociation  of  the  3 p  state  are  not 
unambiguous,  the  lifetime  of  the  state  is  most  likely  limited 
not  by  dissociation  but  by  radiation.  The  predicted  radianve 
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lifetime  of  18  ns  might  well  preclude  the  3 p  state  being  ob¬ 
served  in  any1  ’• 10  of  the  charge  transfer  collisional  reionzia- 
tion  experiments. 

The  energetics  of  neither  the  3s  nor  3 p  Rydberg  states  of 
H50  is  consistent  with  the  observations  of  Gellene  and  Por¬ 
ter1  who  reported  an  ionization  potential  for  D3  O  of  4.3  eV. 
Using  the  experimental  proton  affinity,  from  our  calcula¬ 
tions,  the  lowest  state  is  estimated  to  lie  5.27  eV  below 
H3  0  +  and  the  3p  level  2.92  eV  below  H3  O  + .  Raynor  and 
Herschbach  obtained  values  of  4.68  and  2.80  eV,  respective¬ 
ly,  from  calculated  virtual  orbital  energies  using  Koopman’s 
theorem. 

Not  all  of  the  experimental  observations  on  H3  O  and 
DjO  can  be  completely  explained  by  our  current  knowledge 
of  the  potential  surfaces  of  H30.  Furthermore,  the  correla¬ 
tion  of  the  1 2A '  state  to  the  lowest  OH  +  H2  asymptote  as 
well  as  the  status  of  the  MCSCF-deterrained  minimum  of 
valence  character  on  the  1 2A  *  potential  surface  require 
further  investigation.  Despite  considerable  study  both  ex¬ 
perimentally  and  theoretically,  there  remain  substantial 
gaps  in  our  understanding  of  the  H3  O  system. 
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The  1  'A'.  1  ‘A'.  1 5A',  and  1  ’A"  states  of  LijH  have  been  investigated  at  the  MCSCF/SOCI  level.  The  global  minimum  is  a 
planar  conformation  of  'A'  symmetry.  A  local  minimum  on  the  same  potential  surface  at  a  C)v  pyramidal  geometry,  of  mixed 
ion-pair  and  covalent  character,  is  found  at  a  relative  energy  of  20.30  kcal/mol.  The  barrier  height  for  isomerization  is  predicted 
to  be  1.3  kcal/mol.  The  correlation  diagram  linking  states  of  LijH  to  those  ofLtj+H.  LiH+Li2  and  LijH+Li  is  presented. 


1.  Introduction 

There  is  considerable  intrinsic  interest  in  under¬ 
standing  novel  bonding  mechanisms  in  polyatomic 
molecules.  In  addition,  stable  states  of  molecules  with 
high  energy  content  are  of  potential  practical  signif¬ 
icance  for  the  development  of  new  propellants  [  1  ]. 
Ion-pair  states  of  small  polyatomics  have  been  sug¬ 
gested  as  potentially  energetic  species.  Theoretical 
methods  are  ideally  suited  for  predicting  energies  and 
stabilities  of  modest-sized  molecules.  However,  the 
primary  focus  of  much  theoretical  work  has  been  on 
conventionally  bonded  molecules  and  ion-pair  states 
remain  largely  unexplored. 

We  have  undertaken,  therefore,  a  survey  of  ion-pair 
and  nearby  states  of  polyatomic  molecules  com¬ 
posed  of  light  atoms.  For  the  case  of  H30,  as  re¬ 
ported  in  a  previous  paper  {2],  ion-pair  attraction 
between  Hf  and  O-  does  not  lead  to  a  stable  state. 
The  CJv  symmetry  pyramidal  geometry  with  0"  cen¬ 
tered  over  Hf  in  its  equilateral  triangle  equilibrium 
geometry  is  a  local  minimum  on  the  second  potential 
energy  surface  of  HjO,  but  it  is  unstable  with  respect 
to  dissociation  to  0+H2+H.  A  predicted  [3]  py¬ 
ramidal  local  minimum  on  the  lowest  potential  sur¬ 
face  of  H«  has  also  been  shown  [4,5  ]  to  be  unstable 
with  respect  to  geometric  perturbations.  Both  of  these 
observations  have  been  qualitatively  explained  by 

1  Present  address;  Laboratoire  de  Radioastronomie,  Ecole  Nor- 

male  Superieure,  73005  Paris,  France. 


noting  that  ground  state  neutral  H3  is  not  a  stable 
species.  Therefore,  any  back  charge  transfer  from  0“ 
or  H~  to  Hf  will  lead  to  the  neutral  which  is  re¬ 
pulsive  with  respect  to  H-  +  H. 

In  contrast,  as  suggested  by  Michels  [6],  the  Li3H 
molecule  may  be  more  likely  to  exist  as  a  stable  lon- 
pair  because  both  Lif  and  Li3  are  bound  species. 
Preliminary  calculations  by  Michels  and  Mont- 
gommery  [7  ]  indicate  a  minimum  on  the  lowest  po¬ 
tential  surface  of  Li3H  at  a  C3v  pyramidal  geometry, 
as  expected,  from  the  equilateral  equilibrium  ge¬ 
ometry  of  Lif  for  the  ion-pair  state.  Hydrides  of 
lithium,  in  general,  have  received  substantial  atten¬ 
tion  in  recent  yean  [8-13].  Composed  of  light  at¬ 
oms  and  thus  readily  calculable,  they  can  serve  as 
models  of  metal  compounds  and  of  clusters.  Li3H,  in 
particular,  has  been  included  in  surveys  of  lithium 
hydrides  by  Kato  et  al.  [  8  ]  and  by  Cardelino  et  al. 
[9],  who  predicted  a  stable  species  with  a  planar  ge¬ 
ometry  bound  by  0.58  [8]  (0.48  [9] )  eV  with  re¬ 
spect  to  Li2H+Li,  calculated  at  the  SCF  level. 

In  this  work,  we  have  used  an  MCSCF/CI  ap¬ 
proach  to  investigate  the  lowest  singlet  and  triplet 
states  of  LijH.  Geometries  have  been  optimized  by 
MCSCF  gradient  techniques;  for  each  symmetry  the 
pyramidal  geometries  anticipated  for  ion-pair  states 
have  been  included  in  the  study.  The  correlation  dia¬ 
gram  linking  these  minima  to  different  possible  dis¬ 
sociation  products  has  been  established  and  the  rel¬ 
ative  energetics  provided  within  a  consistent 
computational  model. 
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2.  Calculations 

Because  we  expect  the  ion-pair  states  of  Li3H  to  be 
described  as  Li? -H~,  the  basis  set  for  H  must  be 
sufficiently  diffuse  to  provide  a  good  description  of 
the  anion.  The  4s3pld  Gaussian  basis  set  for  H  con¬ 
sists  of  a  double-zeta  plus  polarization  (DZP)  ba¬ 
sis  *'  augmented  by  two  diffuse  s.(0.07,  0.02),  two 
diffuse  p  (0.09,0.03)  and  one  d  (1.0)  function.  With 
this  basis  set,  the  calculated  electronic  affinity  of  the 
hydrogen  atom  is  0.66  eV,  in  reasonable  agreement 
with  the  experimental  value  of  0.75  eV  (15]. 

For  the  Li  atom,  basis  sets  used  in  previous  work *2 
were  combined  to  give  a  DZP  basis  extended  by  one 
diffuse  s  exponent  denoted  5s2pltL  For  the  ground 
state  of  the  Li  atom,  the  present  basis  yielded  an  SCF 
energy  of  -7.432413  hartree,  compared  to  the  Har- 
tree-Fock  limit  [18],  —7.432726  hartree.  The  adi¬ 
abatic  ionization  potential  of  Li3  calculated  with  this 
basis  set  at  the  second-order  Cl  level  (SOCI)  is  4.06 
eV,  which  may  be  compared  with  the  experimental 
value  of  4.35  ±0.2  eV  [19].  With  the  present 
(5s2pld/4s3pld)  basis,  the  SCF  energy  for  LiH  at 
the  equilibrium  bond  distance  (3.016  Oo)  is 
-7.986199  hartree  as  compared  with  the  near  Har- 
tree-Fock  value  of  Cade  and  Huo  [20]  of 
-7.987313  hartree.  For  Li3H  geometry  optimiza¬ 
tion,  a  double-zeta  Li  basis  was  used 

The  active  space  for  the  MCSCF  calculations  of 
Li3H  included  the  four  orbitals  arising  from  the  Li 
2s  and  H  Is  orbitals.  Analogous  active  spaces  were 
used  for  calculations  on  the  fragment  molecules.  Ge¬ 
ometries  were  optimized  at  the  MCSCF  level  using 
analytic  gradients  and  analytic  second  derivatives. 
For  some  symmetries,  geometries  were  obtained  by 
direct  variation  of  parameters.  Energies  at  the  min¬ 
ima  were  obtained  from  second-order  Cl  (SOCI) 
calculations  that  included  all  single  and  double  ex- 


"  The  two  inner  j  functions  were  taken  from  the  basis  listed  in 
ref.  (14],  The  p  exponent  was  0.9. 

11  The  five  s  functions  were  taken  from  the  first  five  functions 
listed  by  Gerber  and  Schumacher  [16],  The  third  s  was  re¬ 
placed  by  an  exponent  of  0.45.  The  p  and  d  functions  were 
taken  from  the  F,G  basis  set  of  McAdon  and  Goddard  [17], 

"J  For  the  Li  double-zeta  basis  the  two  most  diffuse  s  functions 
were  replaced  by  a  single  function  with  an  exponent  of  0.055 
and  the  d  function  was  omitted. 
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citations  from  the  active  space.  Calculations  were 
carried  out  in  Q  symmetry.  In  ‘A',  ’A",  3A',  and  3A' 
symmetry,  the  SOCI  expansion  resulted  in  16884, 
16274,  22275,  and  22596  configurations,  respec¬ 
tively.  The  GAMESS  [21  ]  and  ALCHEMY  II  [22] 
programs  have  been  used  in  this  work. 


3.  Resalts  and  discussion 

3. 1.  Fragment  geometries  and  energies 

In  order  to  obtain  dissociation  energies  with  re¬ 
spect  to  all  possible  asymptotes  within  a  consistent 
computational  model  and  to  compare  present  work 
with  previous  calculations  and  with  experiment, 
comparable  MCSCF/SOCI  calculations  have  been 
performed  for  the  diatomics  Li2  and  LiH  and  for  the 
triatomics  Li3,  Li2H,  and  Lit . 

Numerous  calculations  of  Li2  and  LiH  have  ap¬ 
peared  in  the  literature;  only  representative  exam¬ 
ples  [13,23]  are  included  in  table  1.  The  discrep¬ 
ancy  of  the  calculated  dissociation  energies  with  the 
experimental  values  is  1.2  and  3.9  keal/mol  for  Li2 
and  LiH,  respectively.  Core  polarization  has  not  been 
included  in  the  present  calculations  [26]. 

The  lowest  two  states  of  Li3,  which  are  predicted 
[9,16,17,27,28]  to  have  C2v  symmetry,  are  nearly 
degenerate  in  energy.  Therefore,  there  has  been  con¬ 
siderable  work  devoted  to  determining  the  ground 
state.  As  shown  in  table  2,  in  agreement  with  pre¬ 
vious  work,  we  predict  a  2B2  ground  state  with  the 
2A|  state  lying  0.56  kcal/mol  higher  in  energy.  Ge¬ 
ometries  for  both  states  are  in  reasonable  agreement 
with  previous  calculations.  The  predicted  atomiza¬ 
tion  energy  of  the  ground  state  lies  3  kcal/mol  out¬ 
side  the  error  ban  of  the  experimental  determina¬ 
tion  [19].  The  optimized  C3,  symmetry  geometry 
for  Li?  is  also  listed  in  table  2. 

The  ground  state  of  the  Li2H  molecule  is  a  2A,  state 
in  C2y  symmetry,  while  the  first  excited  state  is  pre¬ 
dicted  to  be  a  symmetric  linear  species  of  2IU+  sym¬ 
metry,  lying  9.67  kcal/mol  higher  in  energy.  Results 
of  previous  calculations  are  also  summarized  in  ta¬ 
ble  3. 
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Table  1 

Comparison  of  calculated  results  for  Li2  and  LiH  with  experiment  and  previous  calculations 


Uj 

LiH 

present  work 

calc.  (23J 

exp.  (241 

present  work 

calc.  (131 

exp.  (25] 

basis  set 

5s2pld 

12s7p2d 

5s2pld/4s3pld 

4s4p2dlf/4s2pld 

calculation 

MCSCF/SOCI 

CPP+CI*' 

MCSCF/SOCI 

CEPA  b» 

*.<A) 

2.673  c) 

2.676 

2.673 

1.596  *» 

1.599 

1.596 

D,  (kcal/mol) 

23.14 

24.20 

24.36 

54.11 

57.15 

•  57.98 

•’  Core  polarization  potential+CI  (valence).  41  Correlated  electron  pair  approximation.  **  Experimental  value. 

Table  2 

Comparison  of  calculated  results  for  Li? 

.  Lij(2Bj),  and  LiJ(IA, )  states  with  experiment  and  previous  calculations 

Calculations 

Experiment 

[19] 

present  work 

ref.  [27] 

ref.  (16) 

ref.  (28) 

ref.  (17) 

rtf.  (9] 

basis  set 

5s2pld 

4s3p 

6s3pld 

4s2p 

3s2p 

9s5p 

(uncontracted) 

calculation 

MCSCF/ 

valence 

CEPA 

SCF 

GVB/ 

SCF 

soa 

CI+S+D 

+S+D 

CT 

UJTA.) 

Cj,  symmetry 

U 

/  V 

Li  —  Li 

r(A) 

3.099 

3.084 

3.040 

U 

/  Sri 

U  -  U 

C2,  symmetry 

'i 

r,  (A) 

2.820*' 

2.80 

2.768 

2.45 

2.869 

2.851 

ri  (A) 

3.271  •» 

3.18 

3.235 

2.9 

3.352 

3.188 

D,  (kcal/mol)  01 

35.68 

29.41 

34.40 

28.76 

41.5  ±4 

D0  (kcal/mol) 

34.5  «> 

33.25 

4.35  +  2 

IP  (eV)  (adiabatic) 

4.06 

4.0 

3.95 

U,(JAi) 

ClY  symmetry 

/•.(A) 

3.041 

3.07 

3.0  II 

3.138 

3.245 

r:  (A) 

2.756 

2.70 

2.730 

2.761 

2.692 

T.  (kcal/mol) 

0.56 

0 

0.69 

0.19 

2.17 

11  Gradient  optimization.  *"  Dissociation  to  3LL  °  Estimated  from  vibronic  levels  calculated  in  ref.  [  1 6  J. 


3.2.  Li jH  geometries 

Optimized  geometrical  parameters  for  the  1  ‘A’, 
1 3A',  and  1  ■‘A’  symmetries  of  Li3H  are  listed  in  ta¬ 


ble  4  and  energies  are  given  in  table  5.  The  symmetry 
plane  of  C,  symmetry  contains  H  and  Li'  and  bisects 
the  Li— Li  bond,  as  denoted  in  fig.  1 .  The  global  min¬ 
imum  is  a  planar  geometry  of  A'  symmetry,.in  agree- 
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Table  3 

Comparison  of  calculated  results  for  L^H^A, )  and  Li:H(JI^  )  states  with  previous  calculations 


Present  work 

Ref.  [8) 

Rtf.  (9) 

Ref.  [10! 

basis  set 

5s2pld/4s3pld 

S-21G/2s 

9s5p/9s5p 

4s3p/3slp 

calculation 

M'CSCF/ 

SCF 

SCF 

SCF/ 

SCO 

valence  Cl 

JA, 

Cj,  symmetry 

H 

/  \ 

Li  U 

U-H  (A) 

1.730  •> 

1.764 

1.719 

1.720 

Li-Lt  (A) 

2.564 

2.611 

2.514 

2.536 

0,  (kcal/mol)  •’ 

24.22 

18.90 

21.66 

22.4 

U-H-U 

Li-H  (A) 

1.705 

1.684 

1.70 

Te  (kcal/mol) 

9.67 

10.28 

8.50 

*’  Gradient  optimization  Li  basis  5s2p. 

”  With  respect  to  LiH  +  Lt. 

Table  4 

Optimized  geometrical  parameters  °  for  Li,H 

Ref.  (8)  »* 

Ref,  [9]  c> 

1  ‘A" 

Li-Lt 

1.769 

1.781 

1.740 

planar  (C,) 

U-H 

3.963 

4.185 

4.092 

U-U’ 

3.066 

3.239 

3.118 

U-U 

2.617 

2.599 

2.466 

Ref.  [7J 

1  “A1 

Li-U 

2.737 

2.599  (2.664) 

pyramid  (C„) 

U-H 

1.840 

1.886  (1.896) 

d 

0.942 

1.142  (1.109) 

1  ‘A- 

U-H 

1.824 

transition  state 

U’-H 

1.860 

pyramid— planar 

Lt-U’ 

2.708 

U-U 

2.815 

Ref.  [81  61 

1  ’A’ 

Li-H 

1.754 

1.772 

planar  (C,) 

U'-H 

3.561 

3.609 

U-U' 

3.014 

3.085 

U-U 

2.966 

3.028 

1 5A'  •* 

U-H 

1.716 

pyramid  (C,) 

U'-H 

1.858 

U-U' 

3.206 

U-U 

2.549 

*’  Defined  in  fig.  1,  distances  in  A.  bl  Calculation  of  Katoetal.  [8)  described  m  the  third  column  of  table  3. 

°  Calculation  of  Cardelino  et  al.  (9]  descnbed  in  the  fourth  column  of  table  3.  6<3UG*V5CF  (3-21G”/MP2). 
"  One  imaginary  frequency,  277i  cm'1. 
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Table  5 

Calculated  energies  for  Li>H  system 


Total  energy  (hartree) 

Relative  energy al 
(kcal/mol) 

Zero-point  energy  *’ 
(kcal/mol) 

MCSCF 

SOCI 

LijH('A') 

planar 

-22.938538 

-22.974132 

0.0 

4.66 

Li,H(l  'A') 

Cj,  pyramidal 

-22.901177 

-22.941791 

20.30 

3.75 

Li,H(l‘A') 
transition  state 
pyramid— planar 

-  22.90  li  30 

-22.939750 

21.58 

3.72 

LijH(  1 3A') 
planar 

-22.919232 

-22.945625 

17.89 

4  42 

LijH(  1  JA‘ ) 

C,  pyramidal 

-22.905243 

-22.939508 

21.73 

UH+UjIX'S,*)0 

-22.881967 

-22.918979 

34.62 

2.01  d> 

LijH(jA,)+Ljc' 

-22.892861 

-22.921818 

32.83 

LijH(J£r)+Lic' 

-22.882581 

-22.906410 

42.50 

Li,(:B.)+H 

-22.827279 

-22.853793 

75.54 

Ltj(JA,)  +  Hel 

-22.826337 

-22.852909 

76.09 

*’  With  respect  to  LijH(‘A' )  planar  geometry.  SOCI  calculation.  MCSCF calculation  4s2p/4s3p  Id  basis. 
c '  Supermolecule  calculations.  Li2  spectroscopic  value  [  26  J. 


Fig.  I.  Geometry  of  pyramidal  and  planar  conformations  ofLijH. 
The  vertical  distance  from  the  apex  of  the  pyramid  to  the  plane 
of  the  three  Li  atoms  is  denoted  by  d. 


ment  with  previous  work  (8,9).  A  C3»  pyramidal 
structure,  anticipated  for  the  ion-pair  state,  is  found 
to  be  a  local  minimum  with  all  real  vibrational  fre¬ 
quencies,  20.3  kcal/mol  above  the  planar  minimum. 
A  comparison  with  geometries  reported  by  Michels 
and  Montgomery  (7)  for  this  state  is  also  included. 
Differences  up  to  0.2  A  in  bond  length  due  to  dif¬ 
ferences  in  basis  set  and  computational  model  may 
be  noted.  The  transition  state  for  conversion  of  the 
pyramidal  form  to  the  planar  geometry  has  also  been 
determined.  It  is  predicted  to  have  a  low  imaginary 
frequency,  12Gi  cm"1,  and  to  occur  at  a  very  modest 
geometrical  distortion  from  the  pyramid,  with  a  cal¬ 


culated  barrier  height  including  zero-point  energy  of 
1.3  kcal/mol. 

The  lowest  triplet  state,  1 3A\  has  a  planar  min¬ 
imum  geometry  that  lies  17.9  keal/mo!  above  the 
planar  singlet  ground  state.  We  were  not  able  to  de¬ 
termine  a  pyramidal  local  minimum  for  the  3A"  state. 
Restricted  to  C,  symmetry,  the  1 3A  state  has  a  sta¬ 
tionary  point  at  a  pyramidal  geometry  with  one  im¬ 
aginary  frequency,  277i  cm-1.  An  approximate  in¬ 
vestigation  of  the  1  'A”  state  gave  a  pyramidal  (C, 
symmetry)  form  14  kcal/mol  above  the  1 3A’  state. 
The  1 3A’  and  1  ’A"  states  in  planar  geometries  were 
found  not  to  be  stable  with  respect  to  the 
Li2H(2A,  )  +  Li  and  Li,H(2I*  )  +  Li  asymptotes,  re¬ 
spectively,  with  which  they  correlate 

3.3.  Correlation  diagram 

The  correlation  diagram  linking  the  states  ofLijH 
to  those  of  Lij  +  H,  LiH + Lij,  and  Li2H + Li  is  shown 
in  fig.  2.  The  correlations  have  been  drawn  preserv¬ 
ing  the  C,  symmetry  plane  described  above,  con¬ 
taining  the  H  and  Li'  atoms.  Energies  of  the  Li3  +  H 
and  Li2H + Li  asymptotes  and  for  the  lowest  LiH + Li2 
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LljH  L!3K 

Lfg  +  H  UH  <■  U2  Pyramidal  Planar  U2H  +  U 


Fig.  2.  Correlation  diagram  linking  the  sum  of  Li>H  to  those  ofLij  +  H,  LiH  +  Ltj,  and  LijH  +  Li  has  been  drawn  preserving  C,  symmetry 
( see  teat ) .  The  second  state  of  each  symmetry,  shown  by  a  dashed  line,  has  been  calculated  at  the  geometry  of  the  lower  sute. 


asymptote,  listed  in  table  5,  have  been  calculated  for 
the  superraolecule.  Optimized  geometries  for  the 
triatomics,  tables  2  and  3,  have  been  used.  Excita¬ 
tion  energies  for  the  states  of  Li2  were  taken  from 
tabulated  Tc  values  [25]  and  from  the  calculations 
of  Schmidt-Mink  et  al.  [23]  for  the  non-observed 
a3E„ ,  1 3riu,  2  'I^,  and  1  'IT,  states. 

The  lowest  state  of  each  symmetry  is  clearly  bound 
with  respect  to  all.  possible  separated  fragment 
asymptotes.  The  diagram  also  includes  the  second 
state  of  each  symmetry  at  the  geometry  of  the  ground 
state.  The  energies  are  of  qualitative  reliability  only 
because,  except  for  the  2  2A'  state,  these  states  in¬ 
volve  excitation  to  molecular  orbitals  not  included 
in  the  MCSCF  active  space.  They  are,  however,  pre¬ 
dicted  to  be  bound  with  respect  to  the  LiH+Li2 
asymptotes  to  which  they  correlate. 
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4.  Discussion 

One  goal  of  this  work  was  to  investigate  whether 
ion-pair  attraction  in  Li3H,  unlike  the  cases  contain¬ 
ing  Hj"  as  the  positive  ion,  leads  to  a  stable  energetic 
state.  There  is  only  one  state,  a  ‘A’,  arising  from  the 
Lif  ('Ai)+H"  limit,  which  would  lie  at  152  kcal/ 
mol  in  fig.  2  (not  shown).  Analysis  of  the  wave- 
function  of  the  l  ‘  A'  state  at  the  C3»  pyramidal  struc¬ 
ture  indicates  substantial  mixing  of  ion-pair  and  co¬ 
valent  bonding  character.  The  Li-Li  bond  length  of 
2.737  A  is  0.36  A  shorter  than  the  corresponding 
bond  length  in  Lij\  in  contrast  to  the  case  [2]  of 
HjO,  where  the  ion-pair  state  was  found  at  the  H-H 
separation  of  Hj" .  Both  Li-Li  bond  lengths  in  the  C2v 
structure  of  Li3(2A, ),  to  which  this  state  correlates, 
are  longer  as  well. 
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The  prediction  of  a  true  minimum  on  the  1  'A' 
surface  (all  real  frequencies)  at  the  C3v  geometry  of 
mixed  ionic  and  covalent  character  supports  the 
original  suggestion  [6]  that  because  of  the  stability 
of  Lij  and  Lij" ,  back  charge  transfer  will  not  limit  the 
stability  of  Li3H  with  respect  to  dissociation.  Unlike 
the  H4  and  H30  examples  considered  previously  (2- 
5],  there  is  in  this  molecule,  a  lower  energy  planar 
form.  Given  the  very  low  ( 1.3  kcal/mol)  calculated 
barrier  to  interconversion,  the  pyramidal  local  min¬ 
imum  is  predicted  to  be  unstable  relative  to  the  global 
planar  minimum.  While  the  lowest  state  of  Li3H  is 
most  definitely  stable,  it  is  not  particularly  energetic. 
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